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SUSCEPTIBILITY AND RESISTANCE GENES FOR 
BIPOLAR AFFECTIVE DISORDER 



CROSS-REFERENCE TO RELATED APPLICATIONS 

Q This application claims priority to US Provisional Application Ser. No. 

tO 60/062,924, filed October 20, 1997. This application is related to 08/827,568, filed March 
pLO 28, 1997, and 60/014,334, filed March 29, 1996. These disclosures are incorporated herein 
r* by reference for all purposes. 

r BACKGROUND OF THE INVENTION 

U Field of the Invention 

This invention relates to the field of diagnosis and treatment of bipolar 
4 5 affective disorders . 

Background 

The most characteristic features of bipolar affective disorder (manic- 
depressive illness) are episodes of mania (bipolar I, BPI) or hypomania (bipolar II, BPII) that 
are interspersed with periods of depression. If untreated, manic-depressive illness is 
20 associated with an approximately 20% risk of suicide. Even with treatment, this disorder 
constitutes a major public health problem, afflicting approximately one percent of the 
population. Goodwin et al. y Manic-Depressive Illness (Oxford University Press, New York, 
1990). 

Although little is known about the etiology or pathophysiology of bipolar 
25 affective disorder, genetic and environmental factors contribute to its pathogenesis, 
especially in families with multiple affected members. Considerable genetic and 
epidemiologic data derived from twin, family and adoption studies provide compelling 
evidence for a genetic etiology of this disorder, but the mode(s) of inheritance has not been 
identified. Craddock et a/., Ann. Med 25:317-322 (1993). Nonetheless, to date, the majority 
30 of genetic linkage studies of bipolar affective disorder have assumed that it exhibits classical 
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Mendelian inheritance attributable to a single major gene. Segregation analyses have yielded 
inconsistent results with most studies rejecting a single dominant or recessive locus 
inheritance model. However, if only BPI is considered, the best single gene model is 
dominant inheritance. Pauls etaL, Neuropsy. Genet., 60:290-297 (1995), 

Due to the complexities inherent in linkage studies of psychiatric disorders, 
one study has focused on the identification of a gene for bipolar illness in a large Old Order 
Amish pedigree in southeastern Pennsylvania. Egeland et aL 9 Nature, 325:783-787 (1987). 
The Old Order Amish are a religious sect numbering approximately 15,000 who descend 
from some 30 pioneer couples and who have remained genetically isolated, thereby 
minimizing the introduction of multiple genes responsible for inherited disorders. Amish 
families have large sibships and multiple living generations, making them ideal for genetic 
studies. Further, alcohol and drug abuse, which often complicate psychiatric diagnoses, are 
rare among the Amish. Bipolar affective disorder, however, occurs amongst the Old Order 
Amish with a prevalence rate, characteristic symptom pattern and clinical course that are 
similar to those in the general North American population. The identification and 
characterization of these pedigrees led to the initiation of early genetic linkage studies but no 
evidence for linkage between various polymorphic serum proteins or blood group antigen 
loci and affective disorder was found. 

More recently, using a molecular genetic approach, Egeland and colleagues 
reported evidence supporting the localization of a gene conferring a strong predisposition to 
bipolar affective disorder linked to two loci located on the short arm of chromosome 1 1, the 
Harvey-ras-1 oncogene locus (HRAS) and the insulin (INS) locus. Id. However, reanalysis 
of the Old Order Amish pedigree to include several new individuals, two changes in clinical 
status, and a large lateral extension of the original pedigree markedly reduced the probability 
of linkage between bipolar affective disorder and the HRAS and INS loci. Kelsoe et ah, 
Nature, 342:238-243 (1989). 

Attempts to replicate linkage findings for bipolar affective disorder have 
proven problematic and have been plagued by diagnostic uncertainties, genetic 
heterogeneity, phenocopies, genotyping errors, and the complexities of performing and 
interpreting statistical analyses (Egeland et al (1987) Nature 325, 783-787; Pekkarinen et al. 
(1995) Genome Res. 5: 105-1 15; Ginns et al. (1996) Nature Genet. 12, 431-435; NIMH 
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Genetics Initiative Bipolar Group (1997) Am. J. of Med Genetics (Neuropsych. Genetics) 74, 
227-269; Blackwood et al. (1996) Nature Genet. 12, 427-430; Freimer et al (1996) Nature 
Genet. 12, 436-441). Reported linkages of bipolar affective disorder to DNA markers on 
chromosomes 18, 21 and X have been difficult to replicate and several proposed linkages 
5 have been refuted upon reanalysis. Kelsoe et ah, Nature, 342:238-243 (1989), Berrettini et 
al, Proc. Natl Acad Set USA, 91 :5918-5921 (1994), Straub et al, Nature Gen. 8:291-296 
(1994), Baron et al, Nat Genet, 3:49-55 (1993), Pauls et al, Am. J. Hum. Genet, 57:636- 
r - 643 (1995). 

*{,;; Moreover, since the inheritance of BP AD is probably multifactorial, the 

CIO possible involvement of multiple genetic components of small effect and/or the occurrence 
%j of major allelic effects only in epistasis must be considered. In addition to susceptibility 

; s „ ! r 

\~ . alleles, there could be alleles that reduce the risk of developing BPAD in a manner similar to 
* that reported for other complex genetic disorders (Philibert et al. (1997) J. Affective 

basis'. 

5'* Disorders 43, 1-3). If model-based linkage analyses are used, a "false negative" linkage 
[ Ji5 finding could result when individuals inherit disease susceptibility alleles but do not manifest 
: the phenotype due to the presence of "protective" alleles. The inclusion of individuals who 
inherit susceptibility alleles but do not manifest disease because of "protective" alleles, or of 
individuals who inherit "protective" alleles but nevertheless manifest the disease, will also 
reduce the power of model-free (allele-sharing) analyses. Thus, regardless of whether model- 
20 based or model-free analyses are used, "wellness" or "protective" alleles could have a 
significant impact on linkage analyses. 

Given the magnitude of the public health problem associated with bipolar 
illness and the availability of treatments for this disorder, what is needed in the art is a means 
to determine the risk to an individual, who comes from an affected family, of developing 
25 bipolar affective disorder. Given that risk can depend both on susceptibility and protective 
alleles, it is desirable to have means to determine the presence or absence of both types of 
alleles associated with bipolar affective disorder. Quite surprisingly, the present invention 
provides these and other advantages. 
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SUMMARY OF THE INVENTION 

Regions of chromosomes 6, 13, and 15 have been determined to comprise loci 
which are associated with susceptibility to bipolar affective disorder (BPAD), while regions 
of chromosome 4 and chromosome 1 1 are associated with resistance to BPAD. 
5 Compositions and methods to determine the* various forms of these loci are useful for a 
variety of diagnostic procedures. 

In one aspect, the present invention provides genetically based methods and 
_ kits for determining a genotype associated with an increased or decreased susceptibility to 

familial bipolar affective disorder in a family affected by bipolar affective disorder. The 
§;D method comprises determining the genotype of at least one family member, wherein the 

genotype is determined with at least one marker for at least one chromosomal region linked 
:^ to a locus associated with susceptibility to bipolar affective disorder. The chromosomal 
s regions are inclusive of and localized between markers D6S344 and D6S89 on chromosome 
J? 6, markers D13S171 and D13S2 18 on chromosome 13, or markers D15S153 and D15S1 17 
t|5 on chromosome 15, such as at about marker D15S148 on chromosome 15. The bipolar 
CI affective disorder disease status is determined for the family member after the age of onset. 
The genotype and disease status of the family member are compared to determine the 
genotype associated with increased or decreased susceptibility to bipolar affective disorder. 
In one embodiment, the genotype is determined with markers to at least two chromosomal 
20 regions linked to a locus associated with susceptibility to bipolar affective disorder. 
Preferably, the genotype is determined with markers D6S7, D13S1, or D15S45, or 
combinations thereof. In another embodiment, the genotype of an affected family member is 
determined. In a further embodiment, the markers are restriction fragment length 
polymorphisms or microsatellite markers. In yet another embodiment, the genotype which 
25 indicates either the presence or absence of a bipolar illness allele is determined. 

In another aspect, the present invention provides methods and compositions 
for determining the increased or decreased risk of a tested individual developing familial 
bipolar affective disorder by comparing the disease genotype of the tested individual to the 
genotype of a family member which is associated with increased or decreased susceptibility 
30 to bipolar affective disorder. The disease genotype is determined with at least one marker for 
at least one chromosomal region linked to a locus associated with susceptibility to bipolar 
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affective disorder. The chromosomal regions are inclusive of and localized between markers 
D6S344 and D6S89 on chromosome 6, markers D13S171 and D13S218 on chromosome 13, 
or markers D15S153 and D15S117on chromosome 15 such as, for example, at about marker 
D15S148 on chromosome 15. In one embodiment, the genotype of the tested individual is 
5 compared to the genotype of an affected family member. In another embodiment, the 

genotype of the tested individual is determined with markers D6S7, D13S1, or D15S45, or 
combinations thereof. In yet another embodiment, the genotype of the tested individual and 
■ family member are determined at all three chromosomal regions of the present invention, 
tfi The invention also provides genetically based methods and kits for 

hJo determining a genotype associated with an increased or decreased susceptibility to familial 
2l!" bipolar affective disorder in which markers associated with resistance to bipolar affective 
■ ■ 3 disorder are detected. The methods involve determining the genotype of at least one family 
O member, wherein the genotype is determined with at least one marker for at least one 
Z<\ chromosomal region linked to a locus associated with resistance to bipolar affective disorder. 
15 The chromosomal regions are on chromosome 4, inclusive of and localized between markers 
r ,s D4S402 and D4S424 and markers D4S43 1 and D4S404, and on chromosome 11, inclusive 
and localized between D11S394 and D11S29. The bipolar affective disorder disease status is 
determined for the family member after the age of onset. The genotype and disease status of 
the family member are compared to determine the genotype associated with increased or 
20 decreased susceptibility to bipolar affective disorder. In one embodiment, the genotype is 
determined with markers to at least two chromosomal regions linked to a locus associated 
with resistance to bipolar affective disorder. Preferred markers for determining the genotype 
on chromosome 4q include, for example, D4S175, D4S422, D4S1576, D4S2294, D4S1579, 
D4S397, D4S3089, D4S2965, D4S192, D4S420, D4S1644, D4S3334, or combinations 
25 thereof. Preferred markers for determining resistance alleles on chromosome 4p include, for 
example, D4S3007, D4S394, D4S2983, D4S2923, D4S615, AFM a 184za9, D4S2928, 
D4S1065, D4S1582, D4S107, D4S3009, D4S2906, D4S2949, AFM087zg5, D4S2944, 
D4S403, D4S2942 f D4S2984, D4S1602, D4S1511, D4S2311, D4S3048, or combinations 
thereof. On chromosome- 11, preferred markers include, for example, D11S133, D11S147, 
30 CD3D, D11S285, D11S29, or combinations thereof. 



In some embodiments of the invention, the genotype of an affected family 
member is determined. In a further embodiment, the markers are restriction fragment length 
polymorphisms or microsatellite markers. In yet another embodiment, the genotype which 
indicates either the presence or absence of a bipolar illness allele is determined. 
5 In another aspect, the present*invention provides methods and compositions 

for determining the increased or decreased risk of a tested individual developing familial 
bipolar affective disorder by comparing the disease genotype of the tested individual to the 
Q genotype of a family member which is associated with increased or decreased susceptibility 
Zij to bipolar affective disorder. The disease genotype is determined with at least one marker for 
J[t) at least one chromosomal region linked to a locus associated with resistance to bipolar 

affective disorder. The chromosomal regions are on chromosome 4, inclusive of and 
III localized between markers D4S402 and D4S424 and markers D4S431 and D4S404, and on 
f chromosome 11, inclusive and localized between D11S394 and D11S29. In one 
f : embodiment, the genotype of the tested individual is compared to the genotype of an affected 
W5 family member. In yet another embodiment, the genotype of the tested individual and family 
|i member are determined at all three chromosomal regions of the present invention. 

Another embodiment of the invention provides compositions, methods and 
kits for determining the presence of a genotype associated with resistance to bipolar affective 
disorder in a family affected by BPAD. These methods involve determining the genotype of 
20 at least one family member, wherein the genotype is determined with at least one marker for 
at least one chromosomal region linked to a locus associated with resistance to bipolar 
affective disorder. The chromosomal regions are on chromosome 4, inclusive of and 
localized between markers D4S402 and D4S424 and markers D4S431 and D4S404, and on 
chromosome 1 1, inclusive and localized between D11S394 and D11S29. In one embodiment, 
25 the genotype is determined with markers to at least two chromosomal regions linked to a 
locus associated with resistance to bipolar affective disorder. 

In yet another aspect, the invention provides methods and kits for determining 
the contribution of a chromosomal region to the presence or absence of bipolar affective 
disorder, or resistance to BP AD, in a family affected by bipolar affective disorder. The 
30 method comprises determining the corresponding genotype of at least two family members, 
wherein the genotype is determined with at least one marker for at least one tested 
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chromosomal region linked to a locus associated with susceptibility or resistance to bipolar 
affective disorder. The tested chromosomal regions for susceptibility are inclusive of and 
localized between D6S344 and D6S89, D13S171 and D13S218, or at about D15S148; for 
resistance the tested chromosomal regions are inclusive of and localized between either or 
5 D4S402 and D4S424 and markers D4S43 1 and D4S404, and on chromosome 11, inclusive 
and localized between D11S394 and D11S29. The bipolar affective disease status in the 
family members is determined after the age of onset and compared to the genotypes of the 
'-u family members. As a result of this comparison, the contribution of the chromosomal region 
X' = to the presence or absence of bipolar affective disorder in the family is determined. In one 
110 embodiment, corresponding genotype of at least two family members affected by bipolar 
h- illness is determined. In another embodiment, at least one of the markers D6S7, D13S1, or 

D15S45 is used to determine susceptibility, and D4S2949, D4S175, and D4S397 to 
% determine resistance. 

h i BRIEF DESCRIPTION OF THE FIGURES 

T5 Figure 1 summarizes each pedigree in "block" form illustrating that all of the 

BPI pedigrees trace along pathways leading to a common progenitor, one of some 30 couples 
that founded the present Lancaster County, Old Order Amish group. 

Figure 2 shows the maximum lod scores, using BPI as affected diagnosis, for 
two scenarios: (Nl) - nuclear families, homogeneity, dominant inheritance, and affecteds 

20 only; and (N2) - sixteen combinations of analyses for each marker including: dominant vs. 

recessive inheritance, five pedigrees vs. nuclear families, homogeneity vs. heterogeneity, and 
affecteds only vs affecteds and unaffecteds. Nl and N2 represent the number of markers 
furnishing maximum lod scores within given class boundaries for scenarios 1 and 2, 
respectively. 

25 Figure 3 shows the locations of markers on human chromosome 6 that are 

associated with susceptibility to BPAD. The statistical significance of the genetic linkage 
between markers based on sib-pair analysis is shown at left, and map distances between 
markers (in centimorgans) are indicated in the rightmost two columns. 

Figure 4 shows the locations of markers on human chromosome 13 that are 

30 associated with susceptibility to BPAD. The statistical significance of the genetic linkage 



between markers based on sib-pair analysis is shown at left, and map distances between 
markers (in centimorgans) are indicated in the rightmost two columns. 

Figure 5 shows the locations of markers on human chromosome 15 that are 
associated with susceptibility to BPAD. The statistical significance of the genetic linkage 
5 between markers based on sib-pair analysis is shown at left, and map distances between 
markers (in centimorgans) are indicated in the rightmost two columns. 

Figure 6 shows the locations of markers on human chromosome 4p that are 
^ ; associated with resistance to BP AD. The statistical significance of the genetic linkage 
Jt; between markers based on sib-pair analysis is shown at left, and map distances between 
4'0 markers (in centimorgans) are indicated in the rightmost two columns. 
q Figure 7 shows the locations of markers on human chromosome 4q that are 

: associated with resistance to BP AD. The statistical significance of the genetic linkage 
between markers based on sib-pair analysis is shown at left, and map distances between 
markers (in centimorgans) are indicated in the rightmost two columns. 
; 15 Figure 8 shows an analysis of the non-parametric LOD among markers on 

" human chromosome 4q that are associated with resistance to BPAD. 

Figure 9A shows the location of the mouse Clock gene on a genetic map 
chromosome 5 (King et a/., Cell 89: 641-653 (1997)). Figure 9B shows an physical map of 
the mouse chromosome 5 region immediately surrounding Clock. Shown in Figure 9C is the 
20 transcription unit map of the Clock locus. The locations of the homologous region in human, 
which is found on chromosome 4, are indicated. 

Figure 10 shows a summary of the ancestral trace for Amish study bipolar 
pedigrees in "schematic" representation. The LEFT extension coupled with the CORE 
Pedigree 110 provided the resource used to initially report linkage findings (Egeland et al. 
25 (1987) Nature 325, 783-787). Further genetic analyses were reported in 1989 after addition 
of a RIGHT extension to Pedigree 1 10 (Kelsoe et al (1989) Nature 342, 238-243). Pedigree 
210 and partial Pedigree 310 (NIGMS Family 1075) became an additional large lateral 
extension, that along with the earlier subjects, was used in the genome-wide linkage analyses 
reported in 1996 (Ginns et al (1996) Nature Genet 12, 43 1-435). The study reported in 
30 Example IV utilized all of these earlier subjects plus additional expansions, especially in 
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Pedigree 410, so that the overall Study contained 346 samples, including those from 50 BPI 
individuals. 

Figure 1 1 shows a plot of t-statistics obtained from the pair-wise linkage 
results. The figure insert depicts a cumulative plot of p-values whose linearity would reflect 
5 uniformity in /^values associated with multiple linkage results whose null hypotheses were 
all true (see text). The outlying t-statistics and p-values (denoted by arrows) were associated 
with markers, D4S107(t = 6.24), D4S2949 (t = 7.79), D4S2928 (t - 5.03), D11S133 (t - 
6.09), and D11S29 (t = 6 32). 
; : ; Figures 12 A and 12B present a model-free linkage analysis of "wellness" 

do using GENEHUNTER-PLUS -log 10 p. Map position is in Kosambi centimorgans. The -log 10 
- : p was calculated using p values generated by GENEHUNTER-PLUS (including individuals 

- > age 45 yrs in all pedigrees) on the assumption that the NPL score is standard normally 
r f distributed. A -log 10j p of 4.0 corresponds asymptotically to a LOD score of 3.0. Only 
£i: mentally healthy individuals 45 years of age or older were classified as being 'well' (see, 

:f 5 Example IV). Figure 12A: -log 10j p for markers on chromosome 4p: , Pedigree 110 only; 

'9\ and — , Pedigrees 1 10, 210, 3 10 and 410; Figure 12B: -log IOj p for markers on chromosome 
4 q: 9 Pedigree 1 10 only; and , Pedigrees 110, 210, 310 and 410. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

Introduction 

20 In the present invention, regions of chromosome 6, 13, and 15 have been 

identified that comprise loci that are associated with susceptibility to a familial form of 
bipolar affective disorder (BPAD). Ginns et al (1996) Nature Genet 12, 431-435. 
Additional chromosome regions on chromosome 4 and chromosome 1 1 are associated with 
resistance to BP AD. Genotypic identification of the loci associated with either the presence 

25 or absence of familial bipolar affective disorder provides a means to assess the risk of a 

tested individual from an affected family having or developing the disease. Moreover, the 
present invention also provides the means to assess whether these loci are implicated in the 
presence or absence of bipolar illness in an individual. 

Accordingly, the methods and compositions of the present invention provide a 

30 means to alert clinicians to a genetic predisposition towards, or resistance to, bipolar 
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affective disorder. The methods of the invention are useful in genetic counseling of 
individuals from families affected with bipolar illness, and aid in the differential diagnosis of 
bipolar illness from other psychiatric pathologies. 

Definitions 

5 As used herein, "bipolar illness," or "bipolar affective disorder," or "manic 

depression" refer to bipolar I (BPI), bipolar II (BPII), or major depressive disorder (MDD). 
See, "Research Diagnostic Criteria," Spitzer et ah, Arch Gen. Psychtat, 35:773-782 (1978), 
h{ incorporated herein by reference. The term "familial" as applied to the defined terms denotes a 
w genetic contribution to the development of bipolar affective disorder as opposed to a strictly 
jlO environmental etiology. 

7" As used herein "allele associated with increased susceptibility to bipolar illness" 

R? or "bipolar illness allele" or "disease allele" refers to a form of a locus on a chromosome which, 
C ; when present in an individual, directly or indirectly causes or increases the risk of developing 
i J bipolar illness. Similarly, "allele associated with increased resistance to bipolar illness" refers 
| 5 to a form of a locus on a chromosome which, when present in an individual, directly or 
j v indirectly increases the resistance of that individual to bipolar illness. The locus may be any 

DNA sequence, e.g., a gene or genes or fragments thereof or a regulatory element. 

As used herein "locus associated with susceptibility to bipolar illness" refers to 

a locus on a chromosome which in at least one form is an "allele associated with increased 
20 susceptibility to bipolar illness." A "locus associated with resistance to bipolar illness" refers to 

a locus on a chromosome which in at least one form is an "allele associated with increased 

resistance to bipolar illness." 

As used herein, "marker" or "polymorphic marker" refers to a polymorphic 

locus that serves to identify a unique locus on a chromosome. An "informative marker" 
25 appears in different forms on each homologous pair of chromosomes such that inheritance of 

the individual chromosomes can be followed. An informative marker may be comprised of 

two or more markers that individually are not informative. 

As used herein, "family member" refers to an individual's consanguineous 

grandparent, parent, child, or sibling although a more distant blood-relative may be used. The 
30 family member may be alive or deceased. 
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As used herein, "family" refers to two or more consanguineous individuals. A 
family may consist of individuals from the same generation or from 2, 3, 4, 5, 6, 7, 8, 9, or 10- 
15 generations. Thus, a family may consist of ethnic groups or subgroups thereof or a 
geographically secluded interbreeding population having a common ancestor. 
5 As used herein, "linked" refers to the greater association in inheritance of two or 

more non-allelic loci than is to be expected from independent assortment. Loci are linked 
because they reside on the same chromosome. Generally, linked loci are separated by less than 
50 centimorgans, preferably less than 30 or 40 centimorgans, and most preferably less than 20 

Q 

■45 or 10 centimorgans. 

l |) As used herein, "chromosomal region" refers to a length of chromosome which 

tt may be measured by reference to the linear segment of DNA which it comprises. The 5 1 and/or 
r * 3' termini of the chromosomal region can be defined by reference to a unique DNA sequence, 
V i.e., a marker. The chromosomal region may be inclusive or exclusive of the defining 5' or 3* 

terminal DNA sequences. Alternatively or additionally, the 5 f and/or 3' termini of a 
ft 5 chromosomal region can be defined by reference to a length of DNA extending from a unique 
p f DNA sequence. Typically the length extending from a unique DNA sequence is about 10 
^ centimorgans (or million basepairs) or less, and may be 9, 8, 7, 6, 5, 4, 3, 2 or 1 centimorgans 

(or million basepairs) or fractional values thereof wherein the distance in centimorgans is the 

sex-averaged value. 

20 As used herein, "age of onset" refers to the age at which those who develop 

bipolar affective disorder first exhibit its clinically defined symptoms. The age of onset may 
occur at 15 years of age, usually at between 15-20, or 21-25 years of age, and may occur at 26- 
30 or 3 1 -3 5 years of age. 

As used herein, "genotype associated with increased susceptibility to bipolar 
25 affective disorder" refers to a genotype which has a higher probability of occurrence in bipolar 
affective disorder affected family member(s) than in family members who are past the age of 
onset but not affected by bipolar affective disorder. 

As used herein, "genotype associated with increased resistance to bipolar 
affective disorder" refers to a genotype which has a higher probability of occurrence in 
30 individuals who are wholly or partially resistant to BP AD. 

As used herein, a "genotype" may be defined by use of a single or a plurality of 

markers. 
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As used herein, "genotype associated with decreased susceptibility to bipolar 
affective disorder" refers to a genotype which has a lower probability of occurrence in bipolar 
affective disorder affected family member(s) than in family members who are past the age of 
onset but not affected by bipolar affective disorder. 

As used herein, "genotype asspciated with increased or decreased susceptibility 
to bipolar affective disorder" refers to a "genotype associated with increased susceptibility to 
bipolar affective disorder" or a "genotype associated with decreased susceptibility to bipolar 
affective disorder," 

As used herein, "increased" means greater than 50%. 

As used herein, "decreased" means less than 50%. 

As used herein, "determining" the "risk of the tested individual developing 
familial bipolar affective disorder" means ascertaining the probability of the tested individual 
developing bipolar affective disorder after the individual reaches the age of onset. The 
determination of risk may be a quantitatively assessed or may be assessed qualitatively as 
higher, lower, or equivalent to a family member whose corresponding genotype is determined 
at one or more chromosomal regions linked to a locus associated with susceptibility to bipolar 
affective disorder. 

As used herein, "corresponding genotype" refers to a genotype obtained using at 
least one marker from within the same chromosomal region used to genotype another family 
member such that a basis of comparison at that same chromosomal region is provided. A 
corresponding genotype may conveniently be determined using at least one of the same 
markers. 

As used herein, "tested individual" refers to an individual, pre- or post-partum, 
whose genotype is determined and includes a proband. The tested individual is a family 
member from the same family as the family member whose genotype the tested individual's is 
compared to. 

As used herein, "bipolar illness genotype" refers to a genotype determined with 
at least one marker for at least one chromosomal region linked to a locus associated with 
susceptibility to bipolar affective disorder, wherein the tested chromosomal regions are 
inclusive of and localized between D6S344 and D6S89, D13S171 and D13S218, or at about 
D15S148. 
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As used herein, "bipolar illness resistance genotype" refers to a genotype 
determined with at least one marker for at least one chromosomal region linked to a locus 
associated with resistance to bipolar affective disorder, wherein the tested chromosomal regions 
are inclusive of and localized between markers D4S402 and D4S1625 and markers D4S43 1 
5 and D4S404. 

In the form of bipolar affective disorder addressed herein, one or more of the 
loci associated with susceptibility to bipolar affective disorder have a higher probability of 
occurring as a disease allele in a bipolar illness affected family member than in a non-affected 
^ family member. Conversely, in non-affected family members, one or more of the loci which 
4-0 are associated with susceptibility to bipolar affective disorder have an increased probability of 
^; occurring in a form not found in bipolar illness affected family members. This statistical 
>v correlation provides the means of determining whether a particular genotype is associated with 
] increased or decreased susceptibility to bipolar affective disorder. Further, this correlation 
t ? allows one to determine whether and which of the one, two, or three chromosomal regions of 
|1 5 the present invention contribute to bipolar illness in the affected family. And, since 
~f~] susceptibility to bipolar illness increases with the number of bipolar illness alleles of an 

lot* 

-~ h individual, the methods and compositions provide means of determining a tested individual's 

increased or decreased risk of developing bipolar illness. 

Similarly, one or more of the loci associated with resistance to BPAD have a 
20 higher probability of occurring as a resistance allele in a family member that is not affected 

with BPAD than in an affected family member. The statistical correlation provides a means for 

detennining whether a particular genotype is associated with increased or decreased resistance 

to BPAD. 

The methods of the present invention generally comprise determining the 
25 genotype of at least one family member from a family affected by bipolar affective disorder. 
The affected family will have at least one member with bipolar affective disorder, preferably, 
two, three, four, or more members with bipolar affective disorder. As will be clear to those of 
skill in the art, the family affected by bipolar illness will preferably have at least one prior or 
successive generation of family members such that the loci associated with susceptibility to 
30 bipolar illness are transmitted between at least two generations. Accordingly, genotyping of 
two, three, four, or more family members for the bipolar illness genotype is preferred. Even 
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more preferably, these family members will be from two or more different generations; even 
more preferably three or more generations. 

Methods of genotyping are well known to those of skill in the art. Briefly, the 
methods of determining the bipolar illness genotype typically comprise use of at least one 
5 marker for at least one chromosomal region linked to a locus associated with bipolar illness. 
Typically, nucleic acid probes to a marker within these chromosomal region(s) are used for 
genotyping. The markers to the chromosomal regions are sufficiently close to the loci which 
are associated with susceptibility or, depending on the particular chromosomal region tested, 
p resistance to bipolar illness such that following inheritance of the markers allows for following 
j^ft) inheritance of a locus or loci associated with increased or decreased susceptibility or resistance 
'f fl to bipolar affective disorder. Each marker is specific to a chromosomal region and DNA 
C3 sequence variability in markers typically allows a chromosome to be distinguished from its 
jhj i homolog. However, sufficient conservation in DNA sequence by each marker generally allows 
JL„ transmission of the chromosomal region to be traced from generation to generation. A 
C15 statistically significant correlation between the presence or absence of a chromosomal marker 
i s § with the presence or absence of bipolar illness in a family member after the age of onset allows 
H for the determination of the genotype(s) associated with increased or decreased susceptibility or 
resistance to familial bipolar affective disorder. The chromosomal regions of the present 
invention that display linkage to loci associated with susceptibility to bipolar illness are 
20 inclusive of and localized between the markers D6S344 and D6S89 on chromosome 6, 
D13S171 and D13S218 on chromosome 13, or at about D15S148 on chromosome 15, 
generally about 10 centimorgans or 10 million basepairs flaiiking either side of D15S148; 
preferably localized by, and inclusive of at least, marker Dl 5S 1 17. 

Conversely, chromosomal regions of the invention that are linked to loci 
25 associated with increased resistance to BPAD are found on human chromosome 4, more 

particularly on chromosome arm 4p the regions are inclusive of and localized between markers 
D4S43 1 and D4S404 (Figure 6 and Figure 12 A) and on chromosome arm 4q the regions are 
inclusive of and localized between markers D4S402 and D4S1625 (Figure 7 and Figure 12B). 
The chromosomal regions on arm 4p are generally about 10 centimorgans or 10 million base 
30 pairs flanking either side of D4S2949, more preferably about 5 centimorgans flanking either 
side of d4S2949. Examples of suitable markers include, for example, D4S3007, D4S394, 
D4S2983, D4S2923, D4S615, AFM a 184za9, D4S2928, D4S1065, D4S1582, D4S107, 
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D4S3009, D4S2906, D4S2949, AFM087zg5, D4S2944, D4S403, D4S2942, D4S2984, 
D4S1602, D4S1511, D4S2311, D4S3048, or combinations thereof. Particularly preferred 
markers include D4S3009, D4S2906, D4S2949, AFM087zg5, D4S2944, D4S403, D4S2942, 
D4S2984, D4S1602, D4S1511, D4S2311, or combinations thereof. 
5 On arm 4q, the chromosomal^ regions that are linked to loci associated with 

increased resistance to BP AD are typically within about 10 centimorgans on either side of 
D4S397, more preferably within about 5 centimorgans on either side of D4S397. Suitable 
markers include, for example, D4S1 75, D4S422, D4S1576, D4S2294, D4S1579, D4S397, 
2 D4S3089, D4S2965, D4S192, D4S420, D4S1644, D4S3334, or combinations thereof 
j|b An additional chromosomal region that is associated with resistance to BP AD 

^ is found on human chromosome 11. This chromosomal region is inclusive of and localized 
H between markers D11S133 and D11S29. Preferred markers for this region include, for 
Ui example, D11S133, D11S147, CD3D, D11S285, D11S29, or combinations thereof, 
y The genotype or genotypes associated with increased or decreased susceptibility 

115 or resistance to familial bipolar illness is generally determined upon comparison (i.e. 9 
^ correlation) of the genotype of the family member with that family member's bipolar illness 
h " disease status after the age of onset. Comparison of the family member's genotype with the 

family member's disease status allows one to determine the genotype associated with increased 
or decreased susceptibility or resistance to bipolar affective disorder by the use of statistical 
20 methods well known to those of skill in the art. Thus, for example, if the genotype of an 
affected parent and the genotype of an affected child have only one form of an informative 
marker in common, comparison of their disease status with their genotypes implicates the 
particular chromosomal region identified by that common marker as associated with an 
increased risk of developing bipolar illness, or with an increased resistance to genetic and/or 
25 environmentally induced BPAD. Accordingly, the methods of the present invention also allow 
for the formation of pedigrees of sufficient detail such that determination of an allele(s) 
associated with increased susceptibility or resistance to bipolar affective disorder may be 
determined. 

Due to the increased probability of meiotic crossover events between markers of 
30 the present invention and bipolar illness alleles, determining a bipolar illness genotype is 
preferably achieved using closer rather than more distantly related relatives. For similar 
reasons, markers more proximal to the loci associated with increased or decreased 
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susceptibility to bipolar affective disorder are employed, such as D6S7, D13S1, or D15S45, to 
minimize the chance of crossover events. More preferably, two, three, or more additional 
markers flanking D6S7, D13S1, or D15S45 are employed to aid in the detection of a 
recombination event between a marker and the bipolar illness disease allele. Typically, the 
5 markers are separated by 1, 2, 3, 4, or. 5 centimorgans. Preferably, the markers are informative. 

Similarly, for identification of a BPAD-resistant genotype, closer rather than 
more distantly related relatives are preferred, as are markers more proximal to the loci 
associated with increased resistance to BPAD. Such markers on chromosome arm 4p include, 

0 for example, D4S2366, D4S394a, D4S3007, D4S394, D4S2949, D4S1605, D4S1582, 
SO D4S 1 07m, and D4S403 as shown on Fig. 6. On chromosome arm 4q, preferred markers 
J include, for example, D4S422, D4S2423, D4S422a, D4S175, D4S397, D4S3334, and 
H D4S1644 as shown in Fig. 7, On chromosome 1 1, the preferred markers are in the 

fu chromosomal region inclusive of and localized between markers D11S133 and D11S29; 
;;i these include D11S133, D11S147, CD3D, D11S285, D11S29, or combinations thereof. In 
H l 5 each case, the markers are typically separated by 1, 2, 3, 4, or 5 centimorgans. Preferably, the 

1 x I markers are informative. 

The present invention also provides methods and compositions for determining 
a tested individual's increased or decreased risk of inheriting a disease allele. The method 
comprises determining the bipolar illness genotype of a tested individual from the affected 

20 family according to methods described for determining the genotype of a family member. 

Thus, the genotype is determined with at least one marker for at least one chromosomal region 
which is linked to a locus associated with resistance to bipolar illness. The chromosomal 
regions include chromosome arm 4p, where the regions are inclusive of and localized between 
markers D4S431 and D4S404 (Figure 6 and Figure 12 A) and chromosome arm 4q, where the 

25 regions are inclusive of and localized between markers D4S402 and D4S1625 (Figure 7 and 
Figure 12B). An additional region associated with resistance is found on chromosome 1 1 
inclusive of and localized between markers D11S133 and D11S29. Typically, the markers are 
separated by 1, 2, 3, 4, or 5 centimorgans. 

After determining the tested individual's bipolar illness genotype it is compared 

30 to the genotype associated with increased or decreased susceptibility to bipolar affective 
disorder of the affected family. A corresponding genotype is tested such that at least one 
equivalent chromosomal region of the present invention is utilized during comparison of the 
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tested individual's genotype with that of the genotype associated with increased or decreased 
susceptibility to bipolar affective disorder; sometimes two equivalent chromosomal regions are 
compared, often all three chromosomal regions of the tested individual are compared. 
Conveniently, at least one identical marker is used for each equivalent chromosomal region 
5 compared. 

The described comparison provides for a determination of an increased or 
decreased risk of the tested individual developing familial bipolar affective disorder by 
assessing the similarities and differences between the compared genotypes. The absence in the 
€1 tested individual of the form of a susceptibility marker found in the chromosome complements 
1 0 of affected family members signals a reduced risk inheriting a bipolar illness allele and thus, of 
"I!: developing bipolar illness. Conversely, inheritance by the tested individual of a form of the 

susceptibility marker found in affected family members indicates a correspondingly increased 
3 risk of inheriting the bipolar illness allele. Thus, for example, if the same three forms of a 
% : marker are inherited by an affected parent and affected child, the absence of any one of these 
^15 forms of markers in a tested sibling indicates a decreased risk of inheriting the disease allele. In 
p;j contrast, inheritance by the tested sibling of an increasing number of the bipolar illness 
^ genotypes found in the affected family members indicates an increasing risk of inheriting one 
or more disease alleles. A similar analysis applies to testing for increased or decreased risk of 
BPAD because of the absence or presence, respectively, of a chromosomal region that is 
20 associated with an allele that is involved in resistance to BPAD. 

The methods and compositions of the present invention further provide for 
determining whether a chromosomal region of the present invention is, in fact, contributing to 
the presence or absence of familial bipolar affective disorder in a family with at least one 
member affected by bipolar affective disorder. The method comprises determining the 
25 corresponding genotype of at least two family members using methods described for 

determining a genotype associated with increased or decreased susceptibility or resistance to 
familial bipolar affective disorder. Thus, each genotype is determined with at least one marker 
for at least one chromosomal region which is linked to a locus associated with susceptibility or 
resistance to bipolar illness. The chromosomal regions associated with susceptibility are 
30 inclusive of and localized .between D6S34 and D6S89, D13S171 and D13S218, or at about 
D15S148, generally inclusive of a chromosomal region localized by at least D15S1 17. 
Preferably, the markers comprise D6S7, D13S1, or D15S45. More preferably, markers 
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flanking D6S7, D13S1, or D15S45 are also employed. Typically, the markers are separated by 
1, 2, 3, 4, or 5 centimorgans. Chromosomal regions associated with resistance to BPAD are 
generally are inclusive of and localized between D4S402 and D4S424 (Figure 12B); inclusive 
of and localized between D4S431 and D4S404 (Figure 12 A); or inclusive and localized 
5 between D11S394 and D11S29. Preferred markers include, for example, D4S2366, D4S394a, 
D4S3007, D4S394, D4S2949, D4S1605, D4S1582, D4S107m, and D4S403 as shown on Fig. 
6, and D4S422, D4S2423, D4S422a, D4S175, D4S397, D4S3334, and D4S1644 as show in 
Fig. 7. Other preferred markers for resistance include D4S1 75, D4S422, D4S1576, D4S2294, 
§ D4S1579, D4S397, D4S3089, D4S2965, D4S192, D4S420, D4S1644, D4S3334, D4S3007, 
m D4S394, D4S2983, D4S2923, D4S615, AFM a 184za9, D4S2928, D4S1065, D4S1582, 
5 D4S107, D4S3009, D4S2906, D4S2949, AFM087zg5 t D4S2944, D4S403, D4S2942, 
jl ■; D4S2984, D4S1602, D4S1511, D4S23.ll D4S3048, D11S133, D11S147, CD3D, D11S285, 
s and D11S29. The markers are typically separated by 1 , 2, 3, 4, or 5 centimorgans. 

The bipolar affective disorder disease status of the family members may be 
!t5 affected or unaffected, or both. The bipolar affective disease status is assessed for the family 
C members after the age of onset. Corresponding genotypes are determined so that at least one 
marker from within the same chromosomal region is used such that a basis of comparison at 
that chromosomal region is provided. Generally, markers from within two or three different 
chromosomal regions of the present invention are used so that the contribution of these same 
20 chromosomal regions can be determined. Using* statistical methods well known to the skilled 
artisan, the genotypes of the family members are compared to determine if the chromosomal 
region is associated with the presence or absence of familial bipolar affective disorder. A lack 
of a statistically significant correlation between a form of a marker and a particular disease 
status may indicate that the particular chromosomal region identified by that marker does not 
25 contribute to the presence or absence of the disease. The method thereby allows one to exclude 
one, two, or all three chromosomal regions of the present invention from contributing to bipolar 
affective disorder in an affected family. The method may be applied effectively to family 
members from families where bipolar affective disorder is in part genetic, or wholly 
environmental. 

30 The methods of the present invention may be performed on a wide variety of 

human cells including somatic cell hybrids, purified nuclei, chromosomal preparations or 
nucleic acid sequences comprising a marker to a chromosomal region of the present invention. 
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The cells may be somatic or germline and from any time in gestation including fertilized 
embryo or preimplantation blastocysts. Preferably, somatic cells are employed to avoid the 
possibility of meiotic recombination events between a marker and locus associated with 
susceptibility to bipolar illness and to more readily allow determination of the genotype for a 
5 homologous chromosome pair. 

The methods of the present invention may conveniently be practiced with 
informative markers which differ as to sequence or length such as RFLPs (restriction fragment 
length polymorphisms) and microsatellite markers such as STRPs (short tandem repeat 
polymorphisms) or VNTRs (variable number tandem repeats). However, other means to 
^1 0 distinguish between the bipolar illness genotypes may be used, such as but not limited to, 
L; antigenicity, specificity, or activity of encoded proteins or fragments. 
: " " Isolation of nucleic acids from biological samples for use in the present 

rJ invention may be carried out by a variety of means well known in the art. For example, see 
r; those described in Rothbart et al., 1989, in PCR Technology (Erlich ed., Stockton Press, New 
[\S York) and Han et al., 1987, Biochemistry, 26:1617-1625. Kits are also commercially available 
I;? for the extraction of high-molecular weight DNA for PCR. These kits include Genomic 
C, Isolation Kit A.S.A.P. (Boehringer Mannheim, Indianapolis, IN), Genomic DNA Isolation 
System (GIBCO BRL, Gaithersburg, MD), Elu-Quik DNA Purification Kit (Schleicher & 
Schuell, Keene, NH), DNA Extraction Kit (Stratagene, La Jolla, CA), TurboGen Isolation Kit 
20 (Invitrogen, San Diego, CA), and the like. Use of these kits according to the manufacturer's 
instructions is generally acceptable for purification of DNA prior to practicing the methods of 
the present invention. Prior to determining a bipolar illness genotype, the marker or marker 
which defines it may be amplified using such well known amplification means as the 
polymerase chain reaction (PCR) as described in U.S. Patent Nos. 4,683,195; 4,683,202; and 
25 4,965,1 88. In some case, the informative marker may be transcribed into RNA by the cells. In 
this instance, RNA may be used for amplification or for comparison between the tested 
individual and affected family member. 

Of particular use in the present invention as applied to loci associated with 
susceptibility to BP AD are the following: 
30 The primers 5 f -CTCCAGCCTGGGTCACTA-3' (SEQ ID NO:l) and 5'- 

CTAATGC ATGAC AATAATATTTCC A-3 1 (SEQ ID NO:2) which amplify marker D6S344. 
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The clone p7H4 comprising a probe which, with the restriction enzyme EcoRV, 
can define a polymorphism of marker D6S7. Clone p7H4 may be obtained from the American 
Type Culture Collection (ATCC) as purified DNA with the accession number 57429, or as a 
plasmid in E. coli or phage lysate with the accession number 57428. 

The primers 5'-ACCTAAGCGACTGCCTAAAC-3' (SEQ ID NO:3) and 5'- 
CTTGTTCATCTGCCTTGTGC-3' (SEQ ID NO:4) which amplify chromosome marker 
D6S89. 

Also, primers 5-AGTCTCATGTGACAC AAGGCAG-3 ' (SEQ ID NO:5) and 
5'-TGTAACCTGGAAGTAAGGCATG-3' (SEQ ID NO:6) which also amplify marker D6S89. 

The primers 5'-TAGGGCCATCCATTCT-3' (SEQ ID NO:7) and 5*- 
CCTACCATTGACACTCTCAG-3 1 (SEQ ED NO:8) which amplify marker D13S171. 

The clone p7F12 comprising a probe which identifies chromosome marker 
D13S1 . Probe p7F12 is available from the ATCC as purified DNA using accession number 
57007, or in plasmid in E. coli or phage lysate using accession number 57006. Polymorphisms 
can be defined using restriction enzymes Mspl, TaqI, or Bell in conjunction with probe p7F12. 
A region spanning the marker can be amplified with the primers 7F12-Ia 5'- 
TGTAACTATTGGGAGGAAAGA-3' (SEQ ID NO:9) and 7F12-IIa 
5'-TTGTGTAGGACTCTCTAGTTT-3' (SEQ ID NO:10). 

The primers S'-GATTTGAAAATGAGCAGTCC-S* (SEQ ID NO:l 1) and 5'- 
GTCGGGCACTACGTTTATCT-3' (SEQ ID NO: 12) which amplify chromosome marker 
D13S218. 

The probe inserted into clone pEFZ33 which defines an RFLP for chromosome 
marker D15S45 and is available from the ATCC in E. coli or a phage lysate using accession 
number 61006, or as purified DNA using accession number 61007. 

The primers S'-GCACCAACAACTTATCCCAA-S' (SEQ ID NO: 13) and 5'- 
CCCTAAGGGGTCTCTGAAGA-3' (SEQ ID NO:14) which amplify chromosome marker 
D15S117. 

Other probes and primers useful in the present invention are presented in Table 
I. See, e.g., Gyapay et al., "The 1993-1994 Genethon Human Genetic Linkage Map," Nature 
Genet. 7:246-249(1994). . 
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Of particular use in the present invention as applied to loci found on human 
chromosome 4p that are associated with resistance to BP AD are the following: 

The primers 5 ' -AGGC ATACTAGGCCGTATT-3 ' and 5'- 
TTCCC ATC AGCGTCTTC-3 ' , which amplify chromosome markers D4S43 1 and D4S2366; 
5 The primers 5 '-GCTCAC AGAAGTGCCCAATA-3 ' and 5'- 

CCCTGGGTGAAGTTTAATCTC-3 which amplify chromosome marker D4S2935; 
The primers 5 ' -A1TTTTGCTACATTGGTGAC ATA-3 ' and 5'- 
5; CTTC AGGTTCT ACTAGTTC ATGG-3 ' , which amplify chromosome marker D4S3007; 
K The primers 5'-CCCTTGAGCATCCTGACTTC-3' and 5'- 

So G AGTG AGCCCCTGTACTCC A-3 ' , which amplify chromosome marker D4S394; 
H The primers 5 '-ATCAGGGTTCTCCACACAAA-3 ' and 5'- 

?li TTGGTTG AAACTTGTGG AT ATAAA-3 ' , which amplify chromosome marker D4S1 582; 
h The primers 5'-CATTCTAGTAGTTATTGGCTTATCC-3' and 5'- 

tl CAGTTGCTTGATACCTATATTTTTC-3 \ which amplify chromosome marker D4S 1 605; 
; f 5 The primers 5 ' -CCTT ACGGAT AGGGGC AG-3 ' and 5'- 

t'~ CTAATGTCC AGGTCTACGGC-3 ' , which amplify chromosome marker D4S2949; and 

The primers 5 '- AGGTGGCCCTGAGTAGGAGT-3 ' and 5'- 
TTTGAGGGAATGATTTGGGT-3 ' , which amplify chromosome marker D4S403. These and 
additional markers are shown in Table 2. 
20 Of particular use for detecting markers that are associated with resistance to 

BPAD and are found on chromosome arm 4q are the following: 

The primers 5 '-AATGCTTATCTACCAATGAGTG-3 ' and 5'- 
GTGGCTGGGTAGTATTC ATGG-3 ' , which amplify chromosome marker D4S2423; 
The primers 5 '-GGCAAGANTCCGTCTCAA-3 ' and 5'- 
25 TGAAGT AAATTTGGG AG ATTGT-3 ' , which amplify chromosome marker D4S422; 

The primers 5 ' - AGGG AGGTC ATC AGTTC ATT-3 ' and 5'- 
TGTTGC AAACTTTGCTTTTC-3 ' , which amplify chromosome marker D4S397; 

The primers 5 '-TTCTTTGATTCTTCGGGG-3 ' and 5'- 
TTTCTC AGC AAC ATTCCTCT-3 ' , which amplify chromosome marker D4S420; 
30 The primers 5'-TAACATTGACCGCTCCTCTC-3' and 5'- 

C ATCCTTCCTGGTCCCTAGT-3 ' , which amplify chromosome marker D4S1644; 



• 
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The primers 5'-TAAAACTTCTGAATGAAAAG-3 ' and 5'- 
GTAGGGAGGAATAGTTAG-3 ' , which amplify chromosome marker UT2147; 

The primers 5 '-TGC AAACTGTCACTCAAAAG-3 ' and 5'- 
GCCAAGGCTGATCCTC-3 \ which amplify chromosome marker D4S1565; 
5 The primers 5'-GCGCTCTTGGTATATGGTACAG-3' and 5*- 

TGTGGGC AACGTC ACTC-3 ' , which amplify chromosome marker D4S424; and 
The primers 5 ' -GACTCC AAATC AC ATGAGCC-3 ' and 5'- 
r : GTCTCTGC ATTTGCTGGTTT-3 * , which amplify chromosome marker D4S1625. These and 
£;f additional primers that are useful for amplifying chromosomal markers that identify 
ti© chromosomal regions associated with increased resistance to BP AD are shown in Table 3. 

Significantly, the chromosomal region of human chromosome 4q which is 
C : associated with increased resistance to BP AD includes the human homolog of the mouse Clock 
a- gene (Fig. 9 A). Certain alleles of this gene, which is involved in circadian rhythms, are 
t:) implicated by the findings reported herein as being involved in mediating resistance to BPAD. 
if 5 Accordingly, the present invention provides methods of determining a genotype associated 
£i with increased or decreased resistance to familial bipolar affective disorder by determining 
the genotype of an individual using at least one marker for at least one chromosomal region 
linked to the human Clock gene. The chromosomal regions are inclusive of and localized 
between D4S402 and D4S1625. From the genotype, increased or decreased resistance to 
20 bipolar affective disorder is determined. 

The chromosome markers disclosed may be modified by insertions, deletions, 
substitutions, or additions with the proviso that modified sequence be sufficiently 
complementary to identify the same chromosomal markers as the unmodified sequences. As 
will be recognized by those of skill, the complementary sequences of the probes and primers 
25 may likewise be employed or modified. 

The primer pairs for chromosomal markers are also conveniently used as probes 
for the markers. Additional target regions may be identified by walking from known 
chromosome markers as described above. Techniques for chromosome walking are well 
known in the art as described in Sambrook et a/., Molecular Cloning, A Laboratory Manual, 
30 Cold Spring Harbor Press, 1989. Vectors which are optimized for chromosome walking are 
commercially available (e.g., AD ASH and XFIX (Stratagene Cloning Systems, La Jolla, CA). 
New markers may result from physical mapping of the interval defined by markers D6S34 and 
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D6S89, D13S171 and D13S218. Markers D6S7 and D13S1 could serve as convenient focal 
points for mapping of the intervals. Regions proximal to D15S45 may also be vised to identify 
new markers. Those of skill in the art will appreciate that a variety of methods to identify new 
markers may be employed. For example, the chromosomal regions of the present invention 
5 cloned into a yeast artificial chromosome (YAC) library can be identified and isolated by 

identifying the presence of sequences corresponding to the marker sequences identified above. 
Cosmid subclones can be created to provide more detailed physical maps; and AC repetitive 

* '*> hybridization probes could identify additional microsatellite sequences in the cloned regions. 

n ;■■ Other chromosome markers could be used to extend the physical map beyond the boundaries of 

irto the identified markers to yield other markers. 

p Generally, the markers of the present invention will yield directly or indirectly 

J 1 * ( e -g-> upon treatment of a RFLP with a restriction enzyme) at least two distinct bipolar illness 
~». genotypes since one bipolar illness genotype will have been inherited from each parent. In 
W • some cases, however, only one genotype may result if the tested individual received identical 
y .5 forms of the genotype from both parents. In such cases, informative markers providing distinct 
^ genotypes may be used. The sizes of the markers of the tested individual are determined for 
comparison to the size of the markers of the affected family member. Equivalence in size 
between informative markers for the affected family member and tested individual indicates the 
same genotype as defined by that marker. Differences in size between informative markers for 
20 the affected family member indicates different genotypes as defined by that marker. As will be 
understood by the skilled artisan, construction of the pedigree is performed using the methods 
of the present invention to follow the transmission of genotypes associated or not associated 
with bipolar illness as defined by psychological diagnostic criteria. 

Generally, the sizes will be determined by standard gel electrophoresis 
25 techniques as described in Sambrook et ah, Molecular Cloning, A Laboratory Manual, Cold 
Spring Harbor Press, 1989, and Polymeropoulos et ah, Genomics, 12:492-496 (1992). 
Polyacrylaxnide gel electrophoresis is particularly preferred because of its capability of high 
discrimination. Generally, autoradiography is employed to simultaneously visualize and 
identify the markers. Amplification of markers is generally performed with labelled nucleotide 
30 bases that provide a means for identifying the markers following the procedure. Alternatively, 
labelled nucleic acid primers may employed as labelling probes which can hybridize to the 
amplified markers. Typical autoradiographic labels include 32 P, 14 C, 3 H, 125 I, 35 S, or the like. 
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Alternatively, probes may be labelled with visual labels such as photoluminescents, Texas red, 
rhodamine and its derivatives, red leuco dye and 3,3'5,5'-*tetramethylbenzidine (TMB), 
fluorescein and its derivatives, dansyl, umbelliferone and the like or with horse radish 
peroxidase, alkaline phosphatase, or the like. 
5 Although the present invention has been described in some detail by way of 

illustration and example for purposes of clarity of understanding, it will be obvious that certain 
changes and modifications may be practiced within the scope of the appended claims. 

bib Example I 

|G This Example describes the collection of epidemiologic data from pedigree 

IT! members. 

-_ m The genetic-epidemiologic study of bipolar affective disorders among the Old 

ft': Order Amish in southeastern Pennsylvania has been previously described (Egeland et aL, Am. 

? Z, J. Psychiat, 140:56-61 (1983), Egeland, Genetic Studies in Affective Disorders, eds. Papolos & 

IS Lachman, 70-90 (John Wiley & Sons, NY (1994)). Fig. 1 shows that the ancestral line 

encompasses the earliest cases of recurrent, psychiatric illness and the first confirmed cases of 
bipolar affective disorder. Bipolar I disorder among descendants of other pioneers usually 
occurred after intermarriage with the BP progenitor line. (Egeland, Genetic Studies in Affective 
Disorders, eds. Papolos & Lachman, 70-90 (John Wiley & Sons, NY (1994)). On the extreme 

20 left side of the figure, one observes the LEFT extension coupled with the CORE Pedigree 110 
which provided the resource first used to report genetic linkage data (Egeland et al. 9 Nature 
325:783-787(1987)). After follow-up and addition of a RIGHT extension to Pedigree 110, 
further genetic analyses were reported in 1989 (Kelsoe et ah, Nature, 342:238-243 (1989). 
Next, Pedigree 210 and partial Pedigree 310 (NIGMS Family 1075) (Egeland, NIGMS Human 

25 Genetic Mutant Cell Repository, NIH Publication 94-201 1, 408-428, 992-999 (1994)) became 
a second large lateral extension (Pauls et al., Genomics, 1 1 :730-736 (1991)). The present 
report utilizes all of these earlier subjects plus additional expansions, especially in Pedigree 
310. The diagnoses for the 207 individuals in our current linkage study are summarized in 
Table 2. These Old Order Amish kinships continue to provide for lateral and lineal expansion 

30 and have evolved into the IX-Xth generations of descendants at risk. 
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Case ascertainment for mental illness among the Amish began with a 

community- wide network of informants and institutional rosters reviewed with informed 

consent (Hostetter et al., Am. J. Psychiat, 140:62-66 (1983)). Over 400 patient cases have 

been ascertained. A psychiatric review board composed (since 1976) of Drs. James N. Sussex, 

5 Abram M. Hostetter, John J, Schwab, David R. Offord and Jean Endicott used both psychiatric 

interviews (Endicott et aL, Arch Gen. Psychiat, 35:837-844 (1978)) and abstracted medical 

records to perform diagnostic assessments based on strict Research Diagnostic Criteria (RDC) 

C- (Spitzer et aL, Arch Gen Psychiat^ 35:773-782 (1978)). Assessments by this review board 

(i ' ; were made blind to pedigree membership, diagnostic opinions and treatment information in the 

|p medical records, and genetic marker status. As the Board's diagnostic procedures yielded 

confirmed cases of BPI affective disorder, the immediate families of these patients were 

: ; ; evaluated for psychopathology. Pedigree 1 10 was selected (1981) for initial genetic linkage 

L study because of relationships between nuclear families, based on BPI probands, and illness 

7" spanning several generations (Egeland, Genetic Studies in Affective Disorders, eds. Papolos & 

|5 Lachman 70-90, John Wiley & Sons, NY, (1994)). When one examines the relative risk for 

^ individuals used in this linkage study, there is a very high prevalence of affective disorder, with 

age-corrected morbid risk rates for BPI, BPII, and MDD (major depressive disorder) of 17%, 

4%, and 6%, respectively. This gives an overall rate of 27% for major affective disorder in 

these pedigrees. The present sample, which includes extensions to the original family (Fig. 1) 

20 totals 207 members, with 3 1 diagnosed BPI, 50 with other psychiatric diagnoses (Dx), and 126 

unaffected individuals (Table 4). 

TABLE 4: DIAGNOSES FOR THE SAMPLE OF 207 
OLD ORDER AMISH SUBJECTS STUDIED IN GENOME SCAN 

25 PED.110 PED.110 PED.110 PED.210 PED.310 TOTAL 

LeftExt CORE IstRtExt 2ndRt.Ext 

Present 
Sample 



BPI 


3 


11 


4 


2 


11 


31 


BPII 


0 


3 


1 


1 


3 


8 


MDD:recurrent 


1 


5 


3 


2 


4 


15 


MDD:single 


1 


5 


1 


1 


1 


9 


Other Dx, 


0 


9 


1 


3 


5 


18 


AFFECTED 


5 


33 


10 


9 


24 


81 


UNAFFECTED 


5 


52 


21 


19 


29 


126 


GRAND TOTAL 


10 


85 


31 


28 


53 


207 
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Over 125 medical records were abstracted and Board reviewed to document the 3 1 BPI cases. 
The average age of onset for BPI disorder was 22 years. Reliability of the bipolar diagnoses 
5 was checked when 16 of the 31 cases (52%) were evaluated twice, with an average five year 
interval between the blind assessments using different clinical documentation and resulting in 
100% concordance. The high reliability obtained lessens the likelihood of misdiagnoses or a 
false positive BPI in our linkage analyses (Egeland et aU 9 Psychiat Genets 1:5-18 (1990)). 
?fi Apart from RDC diagnoses, the project psychiatric panel also recorded clinical 

t|) opinions in a consensus "clinical diagnosis." There was 100% concordance between these two 
types of diagnostic conclusions (5 board members) for the 3 1 BPI cases and 13 of the 15 cases 
s,, of recurrent major depressive disorder. Of particular interest are the diagnostic results for the 
^ eight cases of BPII. Four were designated BPII by both RDC and clinical opinion. The other 
C four were labelled BPI according to clinical opinion, and two of these actually were classified 
15 as "probable BPI" by the strict Research Diagnostic Criteria. This is important to note because 
^ it shows that true BPII disorder occurs rarely in these pedigrees; BPII appears more as a "BPI" 
f waiting to happen. 

This study of bipolar affective disorder in the Old Order Amish represents a 19 
year longitudinal study of an isolated population in which there is a relatively narrow spectrum 
20 of illness (not one case of schizophrenia occurs in the pedigrees used for linkage analyses) with 
bipolar disorder being the predominant diagnosis. The rigorous longitudinal assessment of 
these Amish pedigrees combined with the systematic and blind psychiatric evaluations and 
diagnoses should also greatly reduce the number of misdiagnoses included in the linkage 
analyses. Moreover, the restricted gene pool characteristic of this relatively closed population 
25 should reduce the number of disease-causing alleles, minimizing the problem of genetic 
heterogeneity. 

Example II 

This Example describes the collection and analysis of genotypic data. 

Genotypic data were collected for 55 1 DNA markers (RFLP and 
30 microsatellite) from 207 pedigree members, including 31 cases of confirmed BPI disorder. 
Blood samples were collected with informed consent and lymphoblastoid cell lines were 
established at the Coriell Institute of Medical Research and/or the National Institute of 
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Mental Health. The NIGMS Human Genetic Mutant Cell Repository catalog contains 
updated pedigree and diagnostic information (Egeland, NIGMS Human Genetic Mutant Cell 
Repository, NIH Publication 94-201 1, 408-428, 992-999 (1994)). DNA was extracted from 
peripheral blood samples and/or immortalized lymphoblastoid cell lines (Neitzel, Hum, 
5 Genet, 73:320-326 (1986)). The RFLP and "microsatellite markers used resulted in a linkage 
map with an average spacing of between 5 and 10 cM (Gyapay et aL, Nature Gen., 7:246- 
249 (1994), Donis-Keller et aL, Cell, 51:319-337 (1987)). Mapping panels were constructed 
to determine the best order of markers typed on the bipolar pedigrees using genotypic data 
3) from the CEPH version 7 database, using the MultiMap linkage analysis program (Matise et 
?fb al, "Automated construction of genetic linkage maps using an expert system (MultiMap): a 
l\ human genome linkage map" Nature Genet. 6:384-390 (1994). Microsatellite markers were 

genotyped individually by previously described methods (Pauls et al., Am. J. Hum. Genet., 
^ 57:636-643 (1995), Sambrook, Molecular Cloning: A Laboratory Manual, 2nd ed., Cold 

Spring Harbor Press, NY (1 989)) and by multiplex procedures adapted from Vignal et al. 
U S "Nonradioactive multiplex procedure for genotyping of microsatellite markers" in Methods 
£ in Molecular Genetics, (ed. Adolph, K.W.) 221-221 (Academic Press, Orlando 1993). In the 
multiplex genotyping procedures, a total of 25 microsatellite markers were analyzed in each 
lane of the genotyping gels. To accomplish this, five markers were coamplified in each PCR 
tube, and five sets of five markers were pooled and precipitated prior to gel loading. Sets of 
20 five microsatellite markers were amplified in 20 yl reactions that included: lOxnM Tris-HCl, 
Ph 8.3, 50 Mm KC1, 0.001% gelatin, 1.5 Mm MgCl 2 , 0.2mM dNTPs, 0.5 units Taq 
polymerase, 1 |aM of each primer (10 primers in total), and 50 ng genomic DNA. Samples 
were denatured at 94°C for 1 minute, followed by 30-35 cycles at 94°C for 15 seconds, 55 C 
for 15 seconds, and 72°C for 15 seconds. After the final cycle, the reactions were incubated 
25 at 72°C for 3 minutes. Following amplification, five sets of PCR amplifications were pooled 
and isopropanol precipitated in 96 well microwell plates. The pooled PCR products were 
resuspended in 10 loading dye containing formamide, denatured at 95°C for 5 minutes and 
loaded onto 6.0% denaturing polyacrylamide gels. After electrophoresis, the gel 
fractionated microsatellite markers were transferred to nylon membranes by capillary 
30 transfer, and visualized by hybridizing membranes with marker specific, chemiluminescent 
probes. One of the oligonucleotides used to amplify each marker was labelled with 
peroxidase using the ECL detection kit (Amersham) and used as probe. Multiple probes 
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corresponding to markers of different sizes were hybridized to the filters simultaneously. 
Chemiluminscent signals were detected by autoradiography. Allele sizes for the 
microsatellite markers were determined relative to a PUC18 sequence or SEQUAMARK 
marker ladder (Research Genetics). To maintain allelic designations for the purposes of 
5 allele frequency calculations, DNA samples "from replicate individuals were included within 
and between gels. Films were scored either manually or using semi-automated allele calling 
software (Biolmage), and were independently analyzed by two individuals blind to disease 

* : : i status. Data from manually scored markers and from the automated scoring system were 

*y transferred into the same file system for linkage analyses. 

f"0 Example III 

This Example describes the statistical analysis of genotypic data. 

Since the exact mode of inheritance of bipolar affective disorder is unknown, 

[ linkage analyses were carried out using nonparametric (allele sharing, model-independent) 

45 methods [SAGE, Sibpal program for Haseman-Elston sib pair tests; affected sibpair analysis 

(ASP) with weighing of multiple affected in the same sibship by number of meioses; affected 

pedigree member (APM); and transmission disequilibrium test (TDT)] 5 as well as lod score 

analyses. 

For nonparametric analyses, based on the asymptotic (theoretical) 
20 distribution of the test statistic, the SIBPAL program furnished formal (asymptotic) p-values 
in the test for an excess proportion of alleles shared IBD (identical by descent) and in the 
Haseman-Elston regression test (Haseman et al., Behav. Genet., 2:3-19 (1972); Keats and 
Elston, Genet Epidemiol Supplement 1:147-152 (1986)). The "true ,f p-value (the empirical 
significance level) is defined as the probability that the observed result or one more extreme 
25 than it is obtained by chance alone. To estimate empirical significance levels associated with 
these results rather than relying on the formal p-value, we carried out computer simulations 
(3000 replicates) for each marker with a formal p-value of 0.01 or less. These simulations 
were extremely time consuming as a complete analysis had to be carried out for each 
replicate. Resulting empirical p-values (identified by * in Table 5) typically were about ten 
30 times higher than the formal "p-values" issued by the program. In addition, a computer 
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program was written to carry out sib pair analyses in which multiple pairs per sibship are 
weighed by the number of meioses. 

Analysis of allele frequencies for the markers D6S7, D13S1 , and D15S45 was 
carried out on probands and their mates to test for the existence of linkage disequilibrium. In 
5 no case was there a significant difference at the 5% level For these three marker loci, as 
well as for the two markers that flanked each, we tested for the presence of specific 
haplotypes. Again, no haplotypes were significantly more frequent than expected by chance 
O when tested at the 5% level. 

fjrj The existence of linkage disequilibrium is known to influence certain types of 

|$ identify-by-state affected sib pair strategies. The basic reason for this is that allele sharing 
P among relatives in isolated populations may be exaggerated due large regions of founder 
i I chromosomes that have not been broken up by recombination. For this reason, we employed 
the TDT to test for linkage, and applied IBD methods, which are unlikely to be significantly 
- = influenced by disequilibrium. 

|| A recent described approach, based on all markers of a chromosome, for 

estimating the proportion of alleles shared IBD was applied to the chromosomes carrying out 
our best three markers. Lander et ah, Nature Genetics, 1 1 :24 1-247 (1995). For allele sharing 
among all possible affected sib pairs, it resulted in suggestive linkage for locus D6S7, The 
TDT did not provide p-values suggestive of linkage when applied to the best three loci and 

20 BPI affecteds. 

For parametric analyses, two-point lod scores were calculated with the 
LINKAGE programs (Lathrop et aU Proc. Natl Acad Set, USA, 81 :3443-3446 (1984)). 
Each marker was analyzed under 16 models (dominant versus recessive inheritance, large 
pedigrees versus data broken into nuclear families, homogeneity versus allowing for 

25 heterogeneity by the HOMOG program (Ott, J., "Analysis of Human Genetic Linkage" 

(Johns Hopkins University Press, Baltimore, 1991)), and for affecteds only versus affected 
and unaffected individuals considered). A penetrance ratio (genetic versus nongenetic cases) 
of 500:1 and disease allele frequencies adjusted to reflect a population prevalence of 1%. 
Individuals without psychiatric illness under any diagnostic scheme were considered 

30 unaffected, whereas those not categorized as affected under one scheme but affected under 
one of the other diagnostic categories were taken to be unknown. Sixteen models were 
tested in the lod score analyses, and included: dominant versus recessive inheritance; large 
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pedigrees versus data broken into nuclear families; homogeneity versus allowing for 
heterogeneity; and analysis for affecteds only versus using affecteds and unaffecteds. To 
evaluate the effect of analyzing the data under multiple models, we compiled maximum lod 
scores for the 551 markers obtained under a single model (case l)(for dominant inheritance, 
5 nuclear families, homogeneity, and affected^ only) and compared them with lod scores 

obtained under multiple testing (case 2) (Fig. 3). Mean lod scores were 0.18 versus 0.46, for 
caseis 1 and 2, respectively. Thus, multiple testing raised lod scores on the average by 0.28 

Q units. Also, no lod score exceeded 2.0 under the fixed model (case 1), whereas 3 lod scores 

m had values of at least 2.0 under multiple testing (case 2). 

Mb Eleven percent of the markers (62 of 551) used in our genome-wide search 

C gave a maximum lod score of 1.0 or higher in at least one of the 16 analysis models. 

fiA Consequently, only regions which yielded stronger evidence of linkage were considered 

lr h further, namely those with markers having a test statistic of p<0.001 in any one analysis type, 

p' 1 or maximum lod score of at least 2.0. Using these criteria, six markers showed evidence for 

|y5 linkage (D1S48, D6S7, D7S67, D11S146, D13S1, and D15S45). 

H Marker Dl 1 S 146 can be obtained from the ATCC using accession Nos. 

59230 for a bacterial/phage lysate, or 5923 1 for purified DNA. Marker Dl S48 (Genethon ID 
G00-000-488) is disclosed in Cartinhour et aL, Cytogenetics and Cell Genetics, 46:591 
(1987). Marker D7S67 (Genethon ID G00-008-432) is disclosed in Donis-Keller et aU Cell 
20 51:319-337(1987). 

Markers at three chromosomal loci gave highly significant test statistic 
p-values under Sibpal (ASP) analyses: D6S7 at chromosome 6pter-p24 with pO.OOOl; 
D13S1 at chromosome 13ql3 with p==0.0003; and D15S45 at chromosome 15qll-qter with 
p=0.0003. Test statistics for these three markers, as well as for markers flanking these 
25 regions, are shown in Table 5. These SIBPAL test statistic p-values were estimated using 
computer simulations (3000 replicates) run under an assumption of no linkage. Therefore, 
they are not flawed by analyses that would furnish spuriously small formal p-values. In lod 
score analysis none of the markers reached LOD=3 criterion. However, in some of the 
nonparametric analysis methods, p-values of less than 0.001 (and even 0.0001, which is 
30 asymptotically equivalent to Z max =3) are found. Together, these results lend further support to 
the significance of these intervals as candidate regions. 
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TABLE 5: RESULTS OF LINKAGE ANALYSES 



Locus 



Map dist Zmax dom Zmax rec SIBPAL p ASP p 



10 



BPI 
D6S344 
D6S70.0 
D6S89 
D6S28 



2.342 
36 . 0 
17.7 



.000 
1.456 
. 097 
.000 



.000 
.0003* 
.001 
.003 



.8729 .2145 
. 0513 

.5567 1.0000 



6262 



7800 



bpi+iii 

2.; D6S344 

t$ D6S7 0 . 0 

u - D6S89 36.0 

D6S28 17.7 



, 000 
.469 
.167 
. 000 



.000 
1.609 
.000 
.000 



9249 
0000 
,7113 
, 6272 



.2341 
.0293 
.3230 
.7867 



r BPI+II 

3© D6S344 

D6S7 0.0 
D6S89 3 6.0 

D6S28 17.7 

2£ BPI+II+MDD 

"1 D6S344 

D6S7 0 . 0 

D6S89 36.0 
D6S28 17.7 



.000 
1.885 
.732 
. 000 



. 000 
1. 606 
.732 
. 000 



,000 
984 
,394 
,000 



000 
795 
399 
000 



9249 
0003* 
5241 
,6874 



8977 
0001 
7126 
6877 



.2036 
,1561 
.2743 
.6892 



.1055 
.4195 
.1714 
.7918 
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TABLE 5: RESULTS OF LINKAGE ANALYSES (cont'd) 
Locus Map dist Zmax dom Zmax rec SIBPAL p ASP p 



10 



25 



30 



35 



40 



45 



BPI 
D13S221 
D13S171 
D13S1 
D13S218 
D13S263 

BPI+II1 
D13S221 
D13S171 
D13S1 
D13S218 
D13S263 

BPI+II 
D13S221 
D13S171 
D13S1 
D13S218 
D13S263 



15 .3 
5.2 
5.1 
10.1 



15 .3 
5.2 
5.1 
10.1 



15 .3 
5.2 
5.1 
10.1 



BPI+II+MDD 
D13S221 

D13S171 15.3 

D13S1 5.2 

D13S218 5.1 

D13S263 10.1 

BPI 
D15S45 

D15S117 5.6 

D15S148 1.2 

D15S38 6.1 

D15S36 0.0 

BPI+II1 
D15S45 

D15S117 5.6 

D15S148 1.2 

D15S38 6.1 

D15S36 0.0 



000 
000 
276 
312 
056 



.000 
.000 
1.402 
.494 
.004 



.000 
.000 
1.203 
.307 
. 006 



.000 
.000 
.000 
.000 
.000 



1.114 
.130 
.338 
.000 
.355 



1.097 
.332 
.752 
.067 
.646 



. 012 
.102 
1.248 
.664 
.175 



.000 
. 000 
1.036 
.423 
.178 



. 000 
. 000 
.676 
.314 
.194 



. 000 
.204 
. 043 
. 025 
. 008 



798 
,580 
,610 
,000 
,400 



,446 
,589 
.613 
,000 
.402 



1. 0000 
.4905 
. 0003* 
.0171 
.0865 



1.0000 
1.0000 
. 0000 
.0175 
.1475 



1.0000 
1.0000 
.0090* 
. 1403 
.2384 



1. 0000 
1.0000 
. 0223 
.3007 
.2262 



. 0003* 
.0843 
.0217 
1.0000 
.0114 



.0018 
.1225 
. 0219 
1.0000 
.0118 



. 9725 
.3736 
.0057 
.0641 
.3028 



.8902 
.4876 
. 0056 
. 0766 
.3998 



.7344 

.4665 

0162 

.1484 

.4033 



.7672 
. 9429 
.0932 
.3202 
.3525 



0163 
.1660 
.1000 
.9853 
. 0862 



0456 
2346 
0976 
9904 
,0844 



50 



BPI+II 
D15S45 
D15S117 
D15S148 
D15S38 
D15S36 



5.6 
1.2 
6.1 
0.0 



.857 
.089 
.461 
.000 
.368 



731 
726 
,829 
, 000 
,292 



.0183 
.0910 
.0123 
1.0000 
. 0131 



0399 
1825 
,0589 
,8551 
,1172 



55 
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TABLE 5: RESULTS OF LINKAGE ANALYSES (cont'd) 



5 


Locus 


Map dist 


Zmax dom 


Zmax rec 


SIBPAL p 


ASP p 




BP+II+MDD 














D15S45 




1.709 4 


.473 


.0032 


.0150 




D15S117 


5.6 


.000 


.096 


.2119 


.3546 


10 


D15S148 


1.2 


.148 


.192 


.0360 


.0998 




D15S38 


6.1 


.000 


. 000 


1.0000 


.7914 




D15S36 


0.0 


.000 


.000 


.1423 


.3168 



Map Dist: Map distance in centimorgan between markers; 
rJ5 Z^: Maximum lod score in analysis of nuclear families, affected only, penetrance ratio (genetic 

Q ; versus nongenetic cases) of 500:1, with allowance for heterogeneity (exception: for D6S7, 

L : „ affecteds and unaffecteds); 

C; ^max^ orn or 2 m ax rec: Under dominant or recessive inheritance; 



SIBPAL p: p-values furnished by SIBAL program in t-test for excess allele sharing in affected sib 

pairs (exception: results for regression analysis given for D6S7). For some markers, an 
.4 empirical p-value, p', was estimated by computer simulation; 

» ASP p: p-values in t-test for excess of allele sharing in affected sibs, multiple sib pairs, in same sibship weighed 
T : by number of meiosis; 

$5 Clinical 

r J Categories: MDD includes only recurrent major depressive 

?\ disorder; Number of affecteds in clinical hierarchies were: 31 BPI, 35 BPI-HBPII, 39 

^ BPI+BPn,and49BPI+BPII+MDD; 

30 1) Only those BP II cases that are borderline BP I are included (such as clinical BP I and RDC manic). 

As observed in Table 5, results are typically stronger for BPI than for more 
liberal diagnostic categories; that is, extending the pool of affected individuals to include 

35 additional psychiatric illness (BPII and recurrent MDD) appears to add "noise" to the 
analyses. Generally, equivalent results are obtained for lod score analyses and our simple 
ASP analysis, whereas the Haseman-Elston approach (SIBPAL program) typically provided 
stronger results. The main differences between the programs SIBPAL and ASP consist in 
the weighing of multiple sib pairs in a sibship (no weighing in SIBPAL). Moreover, 

40 SIBPAL deduces ambiguous genotypes from close or distant relatives while ASP does this 
based only on individuals in the nuclear family. In addition, ASP does not carry out any 
Haseman-Elston type regression analysis as was applied in the case of marker D6S7. The 
fact that some markers flanking our strongly significant markers also show positive linkage 
results provides support for the presence of susceptibility loci near the candidate loci. 
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The relationship between pointwise (locus-specific or nominal) and genome- 
wide significance levels was recently discussed. Lander et al., Nature Genetics, 11:241-247 
(1995)* According to this report, for sib pair methods, pointwise P-values of 0.00074 and 
0.000022 correspond to suggestive and significant linkage, respectively, with "significant" 
5 denoting a genome- wide P-value of 0.05. *For lod score analysis, the respective lod score 
thresholds are L9 and 3.3. Thus, according to these criteria, markers D6S7, D13S1, and 
D15S45 yield locus-specific P-values that are suggestive of linkage. 
p;i Our study of bipolar affective disorder in the Old Order Amish, however, 

Jf represents a 19 year longitudinal study of an isolated population in which there is a 
llO relatively narrow spectrum of illness (not one case of schizophrenia occurs in the pedigrees 
q used for linkage analyses) with bipolar disorder being the predominant diagnosis. The 
k i rigorous longitudinal assessment of these Amish pedigrees combined with the systematic and 
1 blind psychiatric evaluations and diagnoses should also greatly reduce the number of 
-y; misdiagnoses included in our linkage analyses. Moreover, the restricted gene pool 
ri5 characteristic of this relatively closed population should reduce the number of disease- 
Cl causing alleles, minimizing the problem of genetic heterogeneity. 

Similar to other common and complex diseases like diabetes, hypertension 
and perhaps even schizophrenia, our data suggest that genetic factors likely contribute to the 
pathogenesis of bipolar affective disorder, where in the majority of these cases, inheritance is 
20 multifactorial rather than simple Mendelian transmission. Like the genetic variance 
observed for quantitative traits, bipolar affective disorder (even in a relative genetic isolate 
like the Old Order Amish) appears to be a polygenic (complex) trait resulting from the 
variable effects of multiple genes. The results of our genome wide scan suggest that genes 
on chromosomes 6, 13, and 15, rather than just different mutant alleles of a single gene, 
25 determine the susceptibility to and phenotype of bipolar affective disorder in the Old Order 
Amish. Additional sets of genes may underlie the susceptibility to develop bipolar affective 
disorder in other populations. 

Example IV 

This Example describes the ascertainment of psychiatric disorders and health 
30 among several large multigenerational Old Order Amish pedigrees covers a period of over 
twenty years. Throughout this longitudinal study, procedures for assessing and diagnosing 
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subjects have remained constant (Egeland et al. (1990) Psychiat Genet 1, 5-18). Moreover, 
the clinical documentation and diagnostic evaluations have employed rigorous standards and 
been subjected to a variety of reliability tests (Hostetter et al. (1983) Am. J. Psychiat 140, 
62-66). For families in this linkage study, the clinical documentation and diagnostic 
5 evaluations have included a thorough evaluation of all available RDC (Spitzer et al. (1978) 
Arch Gen. Psychiat 35, 773-782) bipolar I (BPI) probands and their relatives. Morbid risk 
analyses have demonstrated a high prevalence of affective disorder among first degree 
% relatives of bipolar probands in these families with the highest risk conferred on the children 
^ of a BPAD parent (Pauls et al (1992) Arch Gen. Psychiat 49, 703-708). Importantly, 
SO because of the long-term, longitudinal nature of the study, even the unaffected, mentally 

healthy individuals (those without any psychiatric illness) in these families have been closely 
1 y followed, many for a period of years past the age of risk for BPAD. Consequently, rather 
p than limit this genome-wide search to identifying susceptibility loci for the disease 
iU phenotype (BPAD), we tested the hypothesis that "protective" alleles may contribute to the 
jL5 absence of psychiatric illness (i.e. mental health "wellness") in unaffected family members in 
H- these "high risk" pedigrees. Since the mode of inheritance of any gene(s) modifying the 

relative risk for affective disorder was unknown (Craddock, N. & McGuffin, P. (1993) Ann. 
Med. 25, 3 1 7-322) we relied exclusively on model-free linkage analyses. 

This Example reports strong evidence for linkage of DNA markers on 
20 chromosome 4p to mental health "wellness" in relatives at high risk for, but who did not 
develop, major affective disorder in several large multigenerational Old Order Amish 
pedigrees with an extremely high incidence of BPAD. 

Materials and Methods 

Diagnostic assessment. 

25 Our genetic-epidemio logic study of BPAD among the Old Order Amish in 

southeastern Pennsylvania has been described in detail (Egeland, J. A. (1994) in Genetic 
Studies in Affective Disorders, eds Papolos, D.F. & Lachman, H.M., (John Wiley & Sons, 
New York) pp. 70-90), including the methods for ascertainment and diagnostic evaluation 
with informed consent (medical records and SADS-L interviews)(Spitzer et al. (1978) Arch. 

30 Gen. Psychiat 35: 773-782; Endicott, J. & Spitzer, R. (1978) Arch Gen. Psychiat 35: 837- 
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844). Diagnoses were made, using strict research diagnostic criteria (RDC)(Spitzer et aL, 
supra.), by a five member psychiatric review board whose members were blind to pedigree 
membership, diagnostic opinions, treatment data from abstracted medical records and genetic 
marker status. By the late 1970's, several dozen BPI probands had been certified by the 
5 psychiatric Board. Subsequently, interviewing began on all available first degree relatives 
using the SADS-L instrument. In this initial screening, over 300 first degree relatives were 
interviewed directly with the SADS-L. These 25 nuclear families, containing one or more 
v cases of BPI, formed the structure of Pedigrees 1 10, 210 and 310 (Figure 10). 
It) The BPI probands in the nuclear families used in this linkage study have on 

lip the average 1 1 .6 first degree relatives. A few siblings were unavailable, while either both 
^" parents (57%) or one parent (23%) were available for interviews and blood samples. Cell 
a lines have been established on an average of eight members for each nuclear family. 

In this study, the "unaffected" individuals (mentally "well" or "healthy") are 
J**] those for whom all SADS-L interview responses were negative (normal) and no 
lt6 contradictory reports were given by family informants. Any individuals for whom some 

symptomatology was identified, even though it did not meet criteria for which the psychiatric 
Board could give a formal diagnosis by RDC, were labeled as "unknowns" in our linkage 
analyses. 

The method used for this longitudinal study is ethnographic and hence 
20 culturally appropriate to the field setting. Each "well" person is not seen annually, nor is 
every individual in a family routinely re-interviewed with the SADS-L. Instead, several 
members of each nuclear family with a BPI proband (BPI nuclear family) are seen annually, 
and those diagnosed with BPI or other major affective disorder undergo a yearly "course-of- 
illness" update. Parents of each BPI patient are regularly visited and they have proven to be 
25 accurate informants about the health of their children and grandchildren. At least one 

"unaffected" sibling (control sample) of the married BPI patients has been interviewed yearly 
since 1990 in connection with a prospective study of "children-at-risk" for bipolar disorder. 
In summary, at least three members and occasionally all members of each BPI nuclear family 
have been evaluated yearly. 
30 Individuals are interviewed anew with the complete SADS-L schedule 

whenever any abnormal mental or emotional symptoms are identified by the follow-up 
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mechanisms. Nearly 50% of those subjects presently carrying a diagnosis of a major 
affective disorder, including BPI, were "unaffected" at the time of the initial SADS-L 
interview. The long-term, systematic follow-up of the families in our study has 
demonstrated that onset of illness in the Old Order Amish is usually reported by multiple 
5 informants. We are confident that individuals designated as "healthy" are free of any 
significant affective disorder. 

v ) Patient samples. 

i?h Blood samples were uniformly collected only after each first degree relative 

Ti - 

12 (including parents, siblings and children older than age 15) of the BPI probands had been 
*i"0 interviewed with the complete SADS-L schedule. Samples were obtained with written 
- U informed consent and coded to maintain confidentiality. The phlebotomist was kept blind to 
^ 3 pedigree relationships and diagnostic status. Lymphoblastoid cell lines were established at 
; 7 ' the Coriell Institute for Medical Research, Camden, N.J. and/or the Clinical Neuroscience 
Ul Branch, IRP, National Institute of Mental Health, Bethesda, MD. The NIGMS Human 
il5 Genetic Mutant Cell Repository catalogue (Egeland, J.A. Amish major affective disorders 
pedigrees. (1994) In 1994A995 Catalog of Cell Lines, NIGMS Human Genetic Mutant Cell 
Repository, 408-428, 992-999 (NIH Publication 94-201 1) contains updated pedigree and 
diagnostic information for several of the Amish pedigrees used in our study. 

Genotyping. 

20 Genomic DNA was obtained from peripheral blood samples and/or 

immortalized lymphoblastoid cell lines as previously described (Ginns et al. (1996) Nature 
Genet 12, 43 1-435). The best order of typed markers on our mapping panels was obtained 
from the genetic location database (LDB) (Collins et al (1996) Proc. Natl. Acad. Sci. USA. 
93, 14771-14775). The order of markers on chromosome 4p is: D4S4l2-6.50cM-D4S431- 

25 0. 24cM-D4S2366-0. 21cM-D4S2935-l. 30cM-D4S3007-1.30cM-D4S394-2. 0cM-D4S2983- 
0. 00cM-D4S2923-0. 00cM-D4S61 5-0. 05cM-AFMal 84za9-L 54cM-D4S2928-L51cM- 
D4S 1065-0. 04cM-D4Sl 582-0. 65cM-D4S107-l, 46cM-D4S3009-0. 30cM-D4S2906-0. OOcM- 
D4S2949-0. 05cM-AFM087zg5-0.24cM~D4S2944-0. llcM-D4S403-0. 40cM-D4S2942- 
0. 00cM-D4S2984-0. 00cM-D4S1602-l.llcM-D4S1511-L 49cM-D4S23ll-2. 15cM-D4S3048- 

30 3.62cM-D4S419~L75cM-D4S404-2.50cM-D4S39L The order of markers on chromosome 
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4q \s:D4S3043-27.91cM-D4S402-0.90cM-D4S427-1.64cM-D4S2303-2.49cM-D4S2985- 
0. 63cM-D4S2423-2. 39cM-D4S2286-l. 50cM~D4S2959-l. 01cM-D4Sl 75-0. 40cM-D4S422- 
0.24cM-D4S1576-4.10cM-D4S2294-0.04cM-D4S1579-0.54cM-D4S397-0.01cM-D4S3089- 
0. 1 0cM-D4S2965-0. 03cM-D4S192-0. 01cM-D4S420-0. 05cM-D4Sl 644-0. 02cM-D4S3334- 

5 0.02cM-D4Sl 565-1 \27chf-D4Sl 625-0 \1 2cM-D4S424-0. Q4cM-D4S 1604-2 SI cM-D4Sl 548. 
The order of markers on chromosome 1 lq is: D11S934-2. lcM-DllS133-8. 7cM-DllS147- 
4. 0cM-CD3D-0. 2cM-Dl 1S285-0. lcM-Dl 1S29. 
% DNA panels for PGR were set up using a 96 microtiter plate format, and the 

PCR master mix was aliquoted using a BioMek robot (Beckman Instruments). PCR was 
HO performed using Perkin-Elmer model 9600 and 9700 thermocyclers. PCR products for a 
M, given DNA marker were optimized by carrying out PCR amplification at 3 different 
\ Kl annealing temperatures on a test panel of genomic DNA samples, and by determining the 

fluorescence signal amplitude and shape following electrophoresis using the ABI 373 
M fluorescent sequencing/genotyping instrument (Applied Biosystems Division, Perkin-Elmer). 
r'l 5 DNA markers were usually processed in groups of six. The genomic DNA samples were 
PCR amplified separately with each of the DNA markers. The PCR products were then 
multiplexed, 6 markers per lane, for electrophoresis on the ABI 373 instruments (Applied 
Biosystems Division, Perkin-Elmer). The DNA from several individuals was represented 
multiple times in the genotyping panels so that within and between each electrophoresis gel 
20 there were "identical" samples that could be used to evaluate the consistency of genotypes 
across several gels. The fluorescent signals from amplified fragments were tracked using 
Genescan (Applied Biosystems Division, Perkin-Elmer), and genotypes were subsequently 
analyzed with Genotyper (Applied Biosystems Division, Perkin-Elmer) . 

Genetic Analysis Software (G.A,S. package version 2.0, Alan Young, Oxford 
25 University, 1 993-1 995) was used to identify problematic marker data, and a utility written in 
SPSS (SPSS Inc.) generated a list of samples that needed to be rerun because of inheritance 
discrepancies or unreadable signals. Samples that had to be rerun were repicked by a 
Microlab 2200 robot (Hamilton Instruments), aliquoted, electrophoresed and analyzed. 
Because we are studying large multigenerational pedigrees where individuals are 
30 descendants of a few progenitors, we maximize the useful information by repeating the 
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genotyping/analysis cycles described above until all possible DNA marker genotypes are 
obtained for the individuals in the study. 

Once genotyping for a marker was finished, the data were reanalyzed with 
G.A.S., observed allelic mutations and other non-inheritances were "zeroed out" in the data 
5 file, and the problematic alleles were notated on pedigree drawings. Histograms were 

generated indicating the marker allele size bins. FASTLINK (Schaffer, A. A. (1996) Hum. 
Hered. 46, 226-235) was used to reanalyze the data prior to further statistical analyses. 

i\ Linkage analyses. 

ft- Model-free linkage analyses were conducted using the two-point affected sib 

rip pair analysis program S.A.G.E. SIBPAL (S.A.G.E. Statistical Analysis for Genetic 

Epidemiology, Release 3.0. (1997) Computer package available from the Department of 
1 Epidemiology and Biostatistics, Rammelkamp Center for Education and Research, 
pi MetroHealth Campus, Case Western Reserve University, Cleveland, OH) and the multipoint 
f; ] analysis program GENEHUNTER-PLUS (Kruglyak, L. & Lander, E.S. (1995) Am. J. 
Jl'5 Hum. Genet 56, 1212-1223). Because there were a few sibships with incomplete marker 
information, marker allele frequencies were estimated from the entire Old Order Amish 
family data set using a maximum likelihood method implemented in the program 
MENDEL AJSERM 1 3 (Lange et al (1988) Genet Epidem. 5, 471-472; Boehnke, M. (1991) 
Am. J. Hum. Genet. 48, 22-25). SIBPAL was used to identify markers showing an excess of 
20 alleles shared identical by descent (IBD) among unaffected, mentally healthy sib pairs. 
Under the null hypothesis of no linkage between a trait and marker, sib pairs would be 
expected to share on the average fifty percent of alleles IBD, but when a trait and marker are 
linked, IBD sharing will be increased in both affected and unaffected sibpairs. Because 
SIBPAL assumes marker allele frequencies appropriate for random samples, it 
25 underestimates the proportion of alleles shared IBD by concordant sib pairs when there is 
linkage. Multipoint analyses using the model-free linkage program GENEHUNTER-PLUS 
produced NPL (non-parametric linkage) scores along points at the chromosomal region of 
interest. Two scoring functions are available in GENEHUNTER-PLUS: IBD sharing can be 
assessed among concordant relative pairs (NPL pairs ) or it may be assessed among larger 
30 groups of concordant relatives (NPL^,). Our analyses were conducted using the wlM statistics 
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as Kruglyak and colleagues have demonstrated that the NPL^ statistic results in a more 
powerful test than the NPL pairs statistic (Boehnke, supra.). 

RESULTS 

First we analyzed our genome- wide scan dataset looking for evidence of 
5 chromosome regions linked to mental health "wellness". In these analyses only mental 
health "wellness" (the absence of any psychiatric illness), in individuals who were over 45 
years of age and had a first degree BPI sibling in their family (Pedigrees 110, 210, 310 and 
2 \ 410), was the linkage phenotype of interest (concordantly unaffected pairs) using SIBPAL. 
5": Of more than 980 DNA markers, only six markers representing three chromosome regions 
JO had t-statistics that were sufficiently outlying and that were likely to represent significant 
t linkage results. Of the markers on chromosome 4p, D4S2949, which is located in the 
vicinity of the BP AD susceptibility locus reported by Blackwood et ah ((1996) Nature 
;; ; ; Genet 12, 427-430), had an empirical SIBPAL p value <5xl 0" 5 (nominal p value <lxl 0~ 7 ). 
H The marker D4S397 on chromosome 4q had an empirical SIBPAL p value = 9x10"* (nominal 
§5 p value = 3xl0' 7 ) On chromosome 1 lq, two DNA markers (D11S133 and D11S29) located 
r over an approximately 20cM region each had a nominal p value <5xl0" 5 (SIBPAL; 
simulations were not performed). To supplement standard criteria for assessing the 
significance of our linkage analysis results, we employed graphical techniques (Figure 11) 
and the empirical assessment of p values (Schweder, T. & Spjotvoll, E. (1982) Biometrika. 
20 69, 493-502; Witte et al (1996) Nature Genet 12, 355-358; Drigalenko, E.L. & Elston, R.C. 
(1997) Genetic Epidem. 14, 779-784). If each marker assessed in a pairwise linkage 
analysis is unlinked to the trait, then the p values associated with those markers should be 
uniformly distributed. In addition, the test-statistics used to generate these p values (for 
instance t-tests in the case of SIBPAL) should follow an appropriate distribution. A plot 
25 (generated using Proc Chart, SAS, SAS Institute Inc.) of the t-statistics obtained from each 
pairwise linkage analysis is shown in Figure 1 1 . The plot in the inset depicts a line that 
should be linear if all markers are unlinked. However, as seen in Figure 11, there are 
outlying t-statistic values that likely represent false null hypotheses; that is, evidence for 
significant linkage results. In addition, in the inset to Figure 1 1, the small upturned portion of 
30 the p value plot near values of l-p = 1 represent departures from uniformity and hence most 
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likely reflect false null hypotheses. Because of the effort required to investigate the 
significance of these findings and the prior evidence supporting a BPAD related locus on 
chromosome 4 (Blackwood et al. (1996) Nature Genet 12, 427-430), we chose to examine 
DNA markers on chromosome 4 first for linkage to mental health "wellness". 
5 To evaluate the findings on chromosome 4p and 4q in more detail, we 

genotyped the subpedigrees and nuclear families containing at least one sibling with BPI 
(Table 6) using additional DNA markers in these interesting regions. Compared to our 
;T;! previous report (Ginns et al. (1996) Nature Genet, 12, 431-435) a larger number of 
^ individuals were included in these analyses (Table 6). In this report, model-free linkage 
Sb analyses using SIBPAL and GENEHUNTER-PLUS (Krugylak et aL (1996) Am. J. Hum. 
Zl Genet 58, 134-1363) were performed using mental health "wellness" as the linkage 
; " phenotype (Tables 7 and 8). In our analyses, individuals having a psychiatric diagnosis other 
P than BPI, as well as those having psychiatric symptoms but no diagnosis, were classified as 
U "unknown category" for affected status. In the Amish Study sample of BPI patients (n = 59) 
}15 the mean and median ages of onset (RDC) are 24 and 22 years, respectively . Hence, in all 
^ analyses we used a conservative age cutoff of 45 years to define family members with the 
unaffected "wellness" phenotype. We also examined the influence of younger age cutoffs 
for defining "well" individuals, and the contribution of different subpedigrees (families from 
pedigrees 1 10, 210, 3 10, and 410 versus only families from pedigree 1 10) on the test 
20 statistics for linkage (Tables 9 and 10). "Well" individuals younger than the specified age 
cutoff were considered to have an "unknown" affected status in the analyses. 

Table 6. Old Order Amish subjects included in linkage analy sis 

Analysis Categories 
Mentally "Unknowns" 
Healthy 

Pedigrees 110, 210, 310, 
410 

k 25 years old 138 85 

^35 years old 109 114 

^45 years old 74 149 
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Analysis 
Mentally 
Healthy 


Categories 
"Unknowns" 


^55 years old 




1/1 


Pedigree 110 only 






^.25 years old 


45 


32 


^35 years old 


37 


40 


2>45 years old 


31 


46 


£.55 years old 


23 


54 



In Table 6, the category of "unknowns" includes individuals of unknown phenotype, 
individuals with psychiatric diagnoses other than BPI, and individuals who are mentally 
healthy but are younger than the particular age cut-off used in analyses. BPI individuals are 
not included in the unknown phenotype category. In pedigrees 110, 210, 310 and 410, 39 
people were diagnosed with BPI, 8 with BPII, 21 with recurrent depressive disorder, 2 with 
unipolar depressive disorder and 15 with other psychiatric illness. In pedigree 110 only, 18 
people were diagnosed with BPI, 2 with BPII, 10 with major depressive disorder, and 5 with 
other psychiatric illness. Note: the individuals used in these linkage analyses represent only 
a subset of the entire Amish bipolar pedigrees since only nuclear families and subpedigrees 
containing a sibling with BPI were included. 

Table 7. Results of SIBPAL analysis of 4p markers 

~ Pedigree 110 Pedigrees 110, 210, 310, 410 

p-value p-value 



Marker 


n(s.e.) 


nominal 


simulated 


ft(s.e.) 


nomi 


simulate 








nal 


d 


D4S412 


.4749(.0621) 


.6555 


np 


.5116(.0539) 


.4154 


np 


D4S431 


.5734(.0441) 


.0523 


np 


.5921(.0388) 


.0110 


np 


D4S2366 


.6781(.0452) 


.0002 


.0005 


.602400356) 


.0027 


.0094 


D4S2935 


.5066(.0218) 


.3825 


np 


.499800198) 


.5043 


np 


D4S3007 


.6233(.0386) 


.0014 


.0023 


.563200337) 


.0330 


.0496 


D4S394 


.6782(.0513) 


.0007 


.0012 


.595500421) 


.0135 


.0249 


D4S2983 


.7219(.0484) 


^xlO" 4 


np 


.609000377) 


.0025 


np 


D4S2923 


.6661 (.0446) 


.0003 


np 


.590200307) 


.0022 


np 


D4S615 


.7161(.0393) 


<lxl0"* 


np 


.622300324) 


.0002 


np 
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Pedigree 110 Pedigrees 110, 210, 310, 410 







P- 


value 




p-value 


Marker 


fl(s.e.) 


nominal 


simulated 


n(s.e.) 


nomi 


simulate 








UAl 


A 
u 


A final 84 


7396( 0446) 


<lxl0' 4 


nn 


6220C 0370^ 


.0008 


np 


xa9 














D4S2928 


.7333(.0257) 


<5xl0" 5 


np 


.636900272) 


<5xl 

n-S 
0 


np 


D4S1605 


.5453(.0258) 


.0440 


.0472 


.579500244) 


.0011 


0.0058 


D4S1582 


6787^ 0616^ 


0032 


.0112 


6269C 0557") 


.0139 


.0510 




6^7 




0090 




<5xl 


.0088 




(.0246) 








O 5 




D4S3009 


.7325(.0552) 


.0001 


np 


.6237(.0379) 


.0008 


np 

XT 




A3Q<t\ 


.vuVrt 


np 


SRSV 0397^ 


.0055 


np 






1 v 1 A-7 


^yI A-5 




<lxl 

o- 7 


<3xl0" 5 


Afin087z 


.522900368) 


.2686 


nD 


.511400246) 


.3218 


np 


D4S2944 


.564700263) 


.0093 


np 


.542800255) 


.0483 


np 


D4S403 


.603200492) 


.0217 


.0233 


.598900443) 


.0232 


.0350 


D4S2942 


.719600308) 


^xlO- 4 


np 


.662700243) 


<lxl 

Q-4 


np 


D4S2984 


.551000396) 


.1032 


np 


.549300297) 


.0505 


np 


D4S1602 


.600100561) 


.0412 


np 


.570300383) 


.0356 


np 


D4S15U 


.624200489) 


.0077 


np 


.577900315) 


.0079 


np 


D4S2311 


.742900279) 


<5xl0" 5 


np 


.632700336) 


.0001 


np 


D4S3048 


.662800573) 


.0036 


np 


.5998(.0403) 


.0078 


np 


D4S419 


.598100270) 


.0004 


.0010 


.577200319) 


.0100 


.0201 


D4S404 


.678500489) 


.0004 


.0010 


.642800470) 


.0020 


.0072 


D4S391 


.7008(.0487) 


.0001 


.0003 


.6585(.0470) 


.0008 


.0035 



II is the estimated proportion of alleles shared identical by descent, 
np: simulations not performed 



Table 8. Results of SIBPAL analysis of 4q markers 



Pedigree 110 Pedigrees 110, 210, 310, 410 

p-value p-value 

Marker f[(s.e.) nominal Simul n(s.e.) nominal simulated 

ated 

D4S3043 .449000531) .8286 np ,5143(.0354) .3442 np 
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Pedigree 110 Pedigrees 110, 210, 310, 410 

p-value p-value 



Marker 


A(s.e.) 


nominal 


Simul 


n(s.e.) 


nominal 


simulate 






ated 






D4S402 


.464900525) 


,/460 


np 


A CnO/ f\A £L'7\ 

.4d9o(.046J) 


QC\A Q 

.o04o 


np 


D4S427 


.475900452) 


H A 1 £L 


«np 


,4j64(.034DJ 


OA/1 A 

.5944 


np 


D4S2303 


.461600423) 




np 


,46/0(.030Dj 


Of Of 
.OJOJ 


np 


D4S2985 


.575400255) 


AAT7 

.002/ 


np 


CMA^/ A1 1A\ 

.MOi^.Oljyj 


A AO C 


np 


D4S2423 


.544500415) 


1 A d 

.1453 


.13 /J 


A CY A1 A/H 

.D44j(.0304j 


.U/4j 


.Uo40 


D4S22S6 


.553300522) 


1 C7A 
.15/0 


np 


.j2ZD(.03oIJ 


.Z/oU 


np 


D4S2959 


.503500359) 


.4619 


np 


.4906(.026oJ 


,63 /U 


np 


T\ AC* 1 TC 

£>4S7 7 J 


rri/A/ AC C 0\ 

.59600055s) 


f\A HI 

.04/1 


.0636 


oyyj(.04o4j 


AOI 1 

.023 1 


A1/1Q 
.UoH-O 


U4IS4J2 


/C1 oo/ a^aan 


A1 HQ 
.0105 


np 




AAA3 


np 


U4oJj/o 




,2oOl 


np 




1 

. 1D4D 


np 


JJ4oJJy4 


.4900^.0440 ) 


.JJJ 1 


n P 




.Ojjo 


np 




A90fV o^ri 


001 S 


tin 

np 




0077 


lip 




7*51 If 0449^ 

. / J 1 1 ^.v l 17 y 


1x1 0" 7 


0009 


6586C 0376^ 


5x1 0" 6 


.0002 




4544^ 0348^ 


.9013 


fit) 


4768f 026 n 


.8120 


np 




^90^ 0^81 




Tin 

"P 


5267C 0366^ 


.2340 


nD 


D4S192 


.5135(.0408) 


.3715 


np 


.504000337) 


.4525 


np 


D4S420 


.5595(.0539) 


.1384 


np 


.5462(.0389) 


.1200 


np 


D4S1644 


.5224(.0521) 


.3351 


.2870 


.5503(.0362) 


.0845 


.0925 


D4S3334 


.5491(.0254) 


.0304 


.0497 


,5258(.0287) 


.1858 


.1769 


D4S1565 


.5091 (.0373) 


.4042 


np 


,504000271) 


.4420 


np 


D4S1625 


.5433(.0454) 


.1730 


np 


.5533(.0339) 


.0603 


np 


D4S424 


.5901(.0527) 


.0481 


np 


.595000461) 


.0226 


np 


D4S1604 


.5501(.0473) 


.1480 


np 


.5095(.0345) 


.3919 


np 


D4S1548 


.5597(.0356) 


.0511 


np 


.5814(.0267) 


.0016 


np 



ft is the estimated proportion of alleles shared identical by descent, 
np: simulations not performed 

Table 9. Results of SIBPAL analysis of selected 4p markers by age 

Pedigree 110 Pedigrees 110, 210, 310, 410 

Marker Age n(s.e.) t-value p-value n(s.e.) t-value p-value 

D4S2366 ^25 .5906(.0415) 2.1813 .0163 .5268(.0220) 1.2195 .1121 
^35 .6295(.0482) 2.6843 .0051 .5459(.0290) 1.5846 .0580 
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Pedigree 110 Pedigrees 110, 210, 310, 410 



Marker 


Age 




t-value 


p-value 




t-value 


p-value 




>45 


.6781 (.0452) 


3.9389 


.0002 


.6024(.0356) 


2.8783 


.0027 


D4S3007 


^25 


.5567(.0341) 


1.6625 


.0505 


.523100196) 


1.1767 


.1205 




^35 


.5767(.0423) 


1.8139 


.0383 


.526300276) 


0.9529 


.1716 




^45 


.6233(.0386) 


3.1978 


.0014 


.563200337) 


1.8729 


.0330 


D4S394 


^25 


.5872(.0397) 


2.1935 


.0159 


.5306(.0220) 


1.3885 


.0834 




*35 


.6242(.0524) 


2.3713 


.0111 


.5560(.0331) 


1.6911 


.0471 




*45 


.6782(.0513) 


3.4734 


.0007 


; 5955(.0421) 


2.2683 


.0135 


D4S1605 


:>25 


.5287(.0259) 


1.1076 


.1361 


.536700227) 


1.6194 


.0541 




^35 


.5271(.0270) 


1.0039 


.1608 


.5405(.0276) 


1.4686 


.0739 




^45 


.5453(.0258) 


1.7551 


.0440 


.579500244) 


3.2623 


,0011 


D4S1582 


k 25 


.5695(.0439) 


1.5849 


.0588 


.507800250) 


0.3110 


.3781 




^35 


.602500586) 


1.7505 


.0436 


.526800358) 


0.7487 


.2280 




^45 


.6787(.0616) 


2.9012 


.0032 


.626900557) 


2.2772 


.0139 


D4S2949 


*25 


.603500305) 


3.3967 


.0006 


.549900205) 


2.4289 


.0081 




*35 


.649700310) 


4.8265 


9x1 0" 6 


.603500260) 


3.9796 


6.8x1 0" 5 




*45 


.707700202) 


10.288 


<lxl0" 7 


.688800243) 


7.7856 


<lxl0" 7 



4 



fl is the estimated proportion of alleles shared identical by descent. 

Table 10. Results of SIBPAL analysis of selected 4q markers by age 



Pedigree 110 Pedigrees 110, 210, 310, 410 



Marker 


Age 


n(s.e.) 


t-value 


P- 


A(s.e.) 


t-value 


p-value 








value 






D4S175 


*25 


.525400402) 


0.6314 


.2650 


.520900283) 


0.7398 


.2305 




;>35 


.5730(0524) 


1.3932 


.0857 


.5557(.0420) 


1,3258 


.0952 




*45 


.5960(.0558) 


1.7200 


.0471 


.5995(.0484) 


2.0536 


.0231 


D4S397 


*25 


.6599(.0307) 


'5.2010 


lxlO- 6 


.630300243) 


5.3595 


2x1 0" 7 




*35 


.7455(.0427) 


5.7485 


6x1 0" 7 


.6622(.0383) 


4.2358 


4.5xl0' 5 




*45 


.751 1(.0449) 


5.5883 


3x1 0- 7 


.6586(.0376) 


4.2156 


5x1 0" 6 


D4S3334 


^25 


.5571(.0265) 


2.1548 


.0174 


.536900187) 


1.9786 


.0246 




5.35 


.5483(.0280) 


1.7269 


.0457 


.5213(.0238) 


0.8968 


.1859 




*45 


.5491 (.0254) 


1.9346 


.0304 


.525800287) 


0.8996 


.1858 
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IT is the estimated proportion of alleles shared identical by descent. 

On chromosome 4p, the maximum multipoint NPL value (GENEHUNTER- 
PLUS) was 4.05 (p = 5.22x1 0" 4 ; including individuals > age 45 yrs in pedigree 110 only) and 
5 4.05 (p = L84xl 0" 4 ; including individuals > age 45 yrs in all pedigrees), respectively. The 
maximum multipoint NPL value (GENEHUNTER-PLUS) for markers on chromosome 4q 
was 3.29 (p = 2.57xl0~ 3 ; including individuals > age 45 yrs in pedigree 1 10 only) and 2.82 (p 
= 4.43x1 0" 3 ; including individuals > age 45 yrs in all pedigrees), respectively. The 
tjl GENEHUNTER-PLUS -log 10j p value as a function of the map position at these locations on 
10 chromosome 4 are shown in Figure 12. SIBPAL test statistics for markers on chromosomes 
£; 4p and 4q are shown in Tables 7 and 8. On chromosome 4 the lowest (nominal) p values 
n[ obtained from the SIBPAL t-statistics were for markers D4S2949 (4p;/?<lxl0* 7 ) and D4S397 
U (4q; p=3xl 0' 7 ). The maximum multipoint NPL value (GENEHUNTER-PLUS) for markers 
1* on chromosome 1 lq was 2.43 (including individuals > age 45 yrs in pedigree 110 only) and 
; 16 2.49 (including individuals > age 45 yrs in all pedigrees), respectively. 
u To obtain empirical /7-values, we simulated genotype data by randomly 

assigning marker alleles to the founders and then assigning alleles to their descendants 
following Mendelian inheritance. Allowing for consanguineous matings, the entire family 
structure (Figure 10) was used in marker assignment, thus taking into account all 
20 relationships between individuals in the dataset. For each simulation, after marker 

assignment, the pedigrees were trimmed down to that of the nuclear families used in the 
linkage analysis. SIBPAL was then run on the trimmed dataset and t-statistics for concordant 
and discordant sib pairs were obtained. The true p value is simply estimated as the proportion 
of replicates in which the simulated statistic is greater than or equal to the observed statistic, 
25 i.e., the probability that the observed result or something more extreme would be obtained by 
chance alone. Simulations were conducted for markers on chromosomes 4p and 4q. For each 
marker, 100,000 replicates were obtained. The empirical p values on chromosome 4p clearly 
meet the proposed criteria of significance for linkage (Lander, E.S. & Kruglyak, L. (1995) 
Nature Genet. 11,241-247). 
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DISCUSSION 

If alleles exist that are associated with mental health "wellness", we reasoned 
that the identification of chromosome regions containing these alleles would be enhanced by 
studying the genetically at risk, mentally healthy members of large, multigenerational 
pedigrees like our Old Order Amish families However, in trying to identify "protective" or 
"wellness" alleles, one must recognize that there are phenocopies that need to be considered. 
Despite the extremely high risk for developing disease, some individuals are undoubtedly 
"well" because they do not inherit any (or all) of the requisite susceptibility alleles for 
BPAD. In addition, since the age of greatest liability for onset of BPAD in the Old Order 
Amish is from early teens through 24 years of age, the misspecification of the "well" 
phenotype for individuals who will eventually develop BPAD would be greatest through this 
age period. In these Old Order Amish families susceptibility alleles for BPAD probably 
occur in very high frequency. Accordingly, an important step in our study which 
demonstrates that there are "protective" alleles was to show that there are "mentally healthy" 
individuals who share marker alleles that should increase the risk of developing BPAD, and 
yet, in the presence of "protective" alleles these individuals do not manifest BPAD. The 
effect of age for inclusion for the "wellness" phenotype can be seen in Tables 9 and 10. For 
many of the markers, ft , an underestimate of the proportion of alleles shared identical by 
descent (IBD) in "well" sibpairs, increases with increasing age, le. a more stringent 
definition of the "well" phenotype. For example, with respect to marker D4S2949 on 4p, ft 
is 0.60, 0.65 and 0.71 for age cutoff points of 25, 35, and 45 years, respectively. This 
suggests that increasing the age for inclusion eliminates some age-related "well" 
phenocopies. 

It is conceivable that virtually all cases of affective disorder in these families 
are due to a common set of susceptibility alleles. The "wellness" or "protective" loci that we 
have tentatively identified could harbor alleles that prevent the manifestation of a bipolar 
affective spectrum disorder phenotype, which could also include major depressive disorder. 
In our analyses the strongest evidence for "protective" alleles comes from pedigree 110, 
suggesting that such alleles may be more likely in this branch of the family. However, 
highly significant test statistics and multipoint lod scores (using GENEHUNTER-PLUS) are 
also observed when pedigrees 1 10, 210, 310 and 410 are used for analyses (Figures 12A and 
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12B). The decreased sharing in proportion of alleles identical by descent (IBD) for 
discordant pairs provides further support for the existence of alleles associated with the 
absence of affective disorder (mental health "wellness") in these families (Table 11). In 
addition, epistatic interactions between alleles could also prevent or delay an illness such as 
5 major depressive disorder from developing into BP AD. Indeed, as we increase the "age of 
risk" cutoff for defining the "well" phenotype from 25 to 45 years in our linkage analyses, 
the number of mentally healthy members decreases as expected, yet the evidence for linkage 
J: increases (Tables 9 and 10). 

rf 0 Table 11* SIBPAL analysis for concordant and discordant pair 



Number of Pedigree 110 Pedigrees 110, 210, 310, 410 

Marker Affected Sibs A P-value A P-value 

(# pairs in 110/alI) ri(s.e.) nominal simulated ft (s.e.) nominal simulated 









CHROMOSOME 4p 










D4S2949 


0 (37/60) 


.7077(.0202) 


<lxl0" 7 


<lx!0" 5 


.6888(.0243) 


<lxl0' 7 


<lxl0" 5 






1 (30/52) 


.5094(.0360) 


.6018 


np 


.4177(.0337) 


.0089 


.0145 






2 (17/20) 


.4183(.0608) 


.9021 


np 


.455900537) 


.7897 


np 


25 






CHROMOSOME 4q 










D4S175 


0 (35/43) 


.596000558) 


.0471 


.0636 


.599500484) 


.0231 


.0348 






1 (27/38) 


.4875(.0611) 


.4194 


np 


.496900513) 


.4762 


np 






2 (17/19) 


.473300528) 


.6901 


np 


.453300533) 


.8042 


np 


30 




















D4S397 


0 (35/43) 


.75il(.0449) 


3xl0" 7 


.0009 


.6586(.0376) 


5x1 0- 6 


.0002 






1 (27/38) 


.4536(.0460) 


.1599 


np 


.506900358) 


.5760 


np 






2 (17/19) 


.500000404) 


.5000 


np 


.5116(.0419) 


.3926 


np 


35 


D4S3334 


0 (37/66) 


.549100254) 


.0304 


.0497 


.525800287) 


.1858 


.1769 






1 (30/56) 


.4119(.0368) 


.0113 


.0089 


.451500317) 


.0655 


np 






2(17/20) 


.445700595) 


.8133 


np 


.4556(.0564) 


.7805 


np 



40 II is the estimated proportion of alleles shared identical by decent. 

Number of affecteds: 0 = mentally healthy (well) sib pairs (older than age 45 years ) 

1 - discordant sib pairs 
2= BPI sib pairs 

np: not performed 



45 



There is some debate on the analysis of sibling pairs as to whether the use of 
inbred sibling pairs results in an increased number of false-positives if allele-sharing-based 
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statistical methods are used (Genin, E. & Clerget-Darpoux, F. (1996) Am. J. Hum. Genet 59, 
1 149-1 162). However, the arguments that a) inbred sibling pairs are likely to share more 
genes than non-inbred sibling pairs (i.e., have a kinship factor greater than 0.5) and b) that 
greater regions of the genome would show significant deviations from the expected non- 
5 inbred sibling sharing value of 0.5, are incorrect when one is merely considering an analysis 
of sibling pairs involving only the transmission of alleles from parents to offspring. The 
transmission of alleles from parents to offspring will follow Mendelian ratios, and thus the 
P , null values for 0, 1, or 2 IBD sibling allele sharing in any population will be 0.25, 0.50, and 
Sf 0.25, whenever only parental and sibling genotype information is used. However, if the 
f§ origin of the parental alleles is taken into consideration, then there will be greater 
p information about alleles shared by sibling pairs from inbred populations. For example, this 

increased informativeness has the potential to resolve ambiguities in the sharing of alleles 
:L transmitted from homozygous parents, since the two copies of the allele in an inbred 
01 homozygous parent could be IBD. This information could also help resolve alleles shared by 
p siblings identical in state into alleles shared IBD, showing that alleles transmitted to two 
F 1 offspring from different parents may be copies of the same allele because of the relatedness 
of the parents. If genealogy is taken into account, then the increased ability to resolve 
ambiguities in allele sharing would result in greater power in the analysis of inbred sibling 
pairs (Genin, E. & Clerget-Darpoux, F., supra.). 
20 Ultimately, if inbreeding exists in a population from which sibling pairs have 

been gathered, but one ignores genealogical information by merely studying the transmission 
of alleles from parents to offspring, then no increase in false-positive linkage results will 
occur. This is because Mendel's law applies to inbred as well as outbred parent-offspring 
allele transmission studies. On the contrary, a decrease in power may result from inbred 
25 sibling pair analyses because spouses may manifest greater homozygosity and therefore 
provide less informative genotypes for parent-offspring-based linkage studies. 

Genetic mapping of complex disorders with multifactorial inheritance could 
be especially difficult if, in addition to susceptibility alleles, individuals inherit "protective" 
alleles that prevent or reduce the risk of manifesting the disease phenotype. Even though 
30 model-based linkage analyses that do not allow for a multifactorial component are of only 
limited usefulness in these circumstances, they are still frequently employed. In these 
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instances, a false negative linkage finding (type 2 error) could result when individuals inherit 
disease susceptibility alleles but do not manifest the phenotype due to the simultaneous 
presence of "protective" alleles. If model-based methods are used, it is important to provide 
a reasonably low estimate of penetrance and include a multifactorial component in the 
model. 

In the initial stages of analyzing a disorder like BPAD which most likely 
displays multifactorial inheritance, robust model-free (allele sharing) methods are usually 
more useful than model-based linkage analysis (Elston, R.C. (1995) Exp. Clin, 
Immunogenet 12, 129-140). Concordant individuals should demonstrate excess allele 
sharing, even with the occurrence of phenocopies, genetic heterogeneity, high frequency of 
susceptibility alleles, and incomplete penetrance. Individuals who inherit susceptibility 
alleles but do not manifest disease because of "protective" alleles, and individuals who 
inherit "protective" alleles but nevertheless manifest the disease will reduce the power of 
these analyses. Thus, regardless of the type of linkage analysis performed, the presence of 
"protective" alleles could have a major impact on identifying susceptibility loci. 

Although the idea that "protective" alleles could modify (or even prevent) a 
behavioral phenotype like BPAD is relatively novel, there are examples where such 
"protective" alleles can affect the expression or inheritance of other Mendelian and 
multifactorial disorders. The severity of sickle cell anemia is influenced by genes that 
increase the amount of circulating fetal hemoglobin (Perrine et aL (1 972) Lancet 2, 1 1 63- 
1 167). Similarly, the genotype of the chemokine receptor CCR5 dramatically influences the 
kinetics of HIV-1 infection, where most individuals who are homozygous for a 32bp deletion 
in the CCR5 gene encoding the coreceptor for macrophage-tropic HIV-1 are "protected" 
from virus infection (Picchio et aL (1997) J. Virology 71, 7124-7127). In Alzheimer^s 
disease, ApoE2, in contrast to ApoE4, appears to reduce the relative risk of developing the 
disease and may protect individuals who inherit a disease-associated ApoE4 allele (Corder 
et aL (1994) Nat Genet 7: 180-183). In an extended Italian family, apolipoprotein A- 
Imilano protects against the development of both clinical and pathologic signs of 
atherosclerosis, despite significantly elevated plasma triglycerides and a markedly decreased 
level of HDL-cholesterol (Franceschini et aL (1980) J. Clin. Invest 66, 892-900). In the non- 
obese diabetic (NOD) mouse model of human autoimmune insulin-dependent diabetes 
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mellitus, partial protection from disease is provided by "resistance" alleles occurring singly 
at either the Idd3 or IddlO non-MHC loci, while epistatic interactions between "resistance" 
alleles at these two loci produces nearly complete protection from diabetes (Wicker et ah 
(1994) J. Exp. Med. 180, 1705-1713). 

There are several mechanisms by which "wellness" or "protective" alleles 
could affect the clinical manifestations of BPAD in the Old Order Amish. One possibility is 
that dominant acting "protective" alleles, either singly or acting together in epistasis, could 
prevent or modify the BPAD phenotype. The variable penetrance of illness or its 
heterogeneous clinical manifestations could result from "resistance" or "protective" alleles 
that alone provide only partial protection, while together with other genes produce epistatic 
interactions resulting in a greater degree of modification of the phenotype. Alternatively, 
there also could be cellular target molecules, e.g. mood "effectors", having forms that are 
either resistant or susceptible to the genetic and/or environmental susceptibility factors for 
BPAD. Individuals having "resistant" mood effectors would be protected from the effects of 
susceptibility alleles and/or environmental factors that result in the BPAD phenotype. In 
contrast, individuals with "sensitive" forms of these mood effectors would be vulnerable to 
developing the BPI phenotype when requisite BPAD susceptibility alleles and/or 
environmental factors are present. 

If epistatic interactions are required for manifestation of the effects of either 
susceptibility or "protective" alleles, the existence of "resistant" and "sensitive" forms of 
cellular effectors or "protective" alleles would be most apparent in families (or populations) 
where there is a high density of affected individuals such as the Old Order Amish in the 
present study. Regardless of the mechanism, the presence of "wellness" or "protective" 
alleles can have a significant impact on linkage analyses as evidenced by preventing the 
appearance of the BPAD phenotype (or its presentation as a forme fruste) in individuals who 
are otherwise genetically predisposed to developing illness. 

Accordingly, a multilocus approach that considers both additive and 
subtractive influences of alleles on the BPAD phenotype is preferred in the identification of 
chromosomal loci harboring genes that contribute to the clinical manifestations of BPAD. 
The involvement of "protective" or "wellness" alleles in determining the manifestation of the 
BPAD phenotype provides an attractive explanation for at least some of the difficulty 
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encountered in searches for BPAD susceptibility alleles. The test statistics from our analyses 
for alleles linked to the absence of psychiatric illness in the Old Order Amish are at least as 
significant as those reported for any susceptibility locus. The identification and 
characterization of "protective" alleles and their gene products can lead to the development 
of a more rational and direct approach to effective therapy for affective disorders. 

All publications and patents mentioned in this specification are incorporated 
herein by reference into the specification to the same extent as if each individual publication 
or patent was specifically and individually indicated to be incorporated herein by reference. 



